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I .  INTRODUCTION 


The  bonding  of  polymers  to  metals  is  accomplished  through  a 
combination  of  chemical  bonding  and  physical  interlocking  of  the  polymer 
around  surface  asperities  on  the  metal.  Many  detailed  studies  of  this 
process  have  been  conducted  {see  Refs.  1-3  and  the  citations  therein).  It 
is  well  known  that  the  metal  interface  must  remain  in  an  uncorroded  state 
in  order  to  preserve  adhesion  during  component  aging.  Much  less  is  known 
about  how  to  preserve  that  metal /polymer  interface  in  the  presence  of 
neighboring  material  layers  that  evolve  liquid  or  gaseous  by-products 
during  aging.  Examples  are  (1)  thermal  insulators  that  must  be  held  on 
spacecraft  outer  bond  lines  by  siloxanes  that  evolve  acetic  acid  during 
cure;  (2)  the  various  layers  of  coatings,  liners,  insulators,  and 
propellants  that  must  adhere  to  solid-rocket  motor  cases  despite  the  high 
reactivity  of  some  of  the  constituent  materials;  and  (3)  the  adhesive 
"flypaper"  used  in  some  electronics  packages  in  order  to  catch  stray  solder 
particles  and  hold  them  for  long  periods  of  time. 

It  is  desired  to  investigate  the  effect  of  filled  polymer  liners  in 
retarding  diffusion  of  corrosive  species  to  the  metal/polymer  interface. 
Consider  the  etched  and  bondable  metal  surface  shown  in  Fig.  1-1.  After 
corrosive  attack  [this  report  concerns  gaseous  hydrogen  fluoride  (HF)  as  an 
example  of  a  corrosive  species],  the  same  metal  surface  has  the  appearance 
shown  in  Fig.  1-2.  Clearly,  a  barrier  liner  on  this  metal  that  could 
prevent  interaction  with  the  HF  would  be  desirable.  The  general  situation 
considered  in  this  report  is  shown  in  Fig.  1-3. 

The  transport  coefficients  of  proposed  liner  materials  must  be 
determined  before  calculations  can  reveal  the  extent  to  which  the  barrier 
reduces  the  flux  (J2  in  Fig.  1-3)  of  corrosive  species  (here,  HF  gas)  at 
the  metal  substrate.  Therefore,  in  Section  II  it  will  be  shown  how  these 
transport  coefficients  are  determined  f  1  om  one-dimensional  diffusion 
experiments.  Concentration  profiles  are  measured  in  order  to  determine  the 
basic  diffusion  coefficient,  D2,  and  a  second  parameter,  H,  which  specified 
the  degree  of  reactivity  between  the  gas  and  the  liner  material.  In 
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Fig.  1-3.  Definition  of  Geometry  and  Material  Parameters.  The 

concentrations,  Cj^x,  t),  of  corrosive  particles  give  rise  to 
fluxes,  J^,  of  corrosive  particles  impinging  on  the  metal 
substrate  and  on  the  inert  boundary  material.  The  matrix 
material  from  x  =  -l1  to  x  =  0  has  a  diffusion  coefficient,  D1f 
and  evolves  corrosive  particles  at  a  rate  r  =  Ge”8t.  The 
barrier  liner  from  x  =  0  to  x  =  l1  has  a  diffusion  coefficient, 
D2,  and  reacts  with  the  corrosive  species  with  a  characteristic 
rate,  H. 
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addition,  the  diffusion  coefficient,  Dj,  of  the  corrosive  species  in  the 
original  matrix  must  be  known  or  determined,  as  well  as  the  rate  of 
corrosive- species  evolution,  r.  We  take  r  =  G  exp(-gt),  where  G  and  g  are 
constants  and  t  is  time. 

In  Section  III,  several  sets  of  diffusion  coefficients  and 
reactivities  are  assumed  for  the  two-layer  system  shown  in  Fig.  1-3.  Time- 
dependent  concentration  profiles  and  cumulative  fluxes  are  calculated.  A 
barrier  effectiveness  is  defined  which  is  zero  when  half  the  number  of 
evolved  particles  impinges  on  the  metal  substrate  (the  other  half 
disappearing  in  the  sink  material  on  the  side  of  the  decomposing  layer 
opposite  to  the  metal  substrate).  The  barrier  effectiveness  is  100?  when 
none  of  the  evolved  particles  impinges  on  the  metal  substrate.  For  the 
case  of  an  0.277 -cm  thick  decomposing  layer  next  to  an  0.025-cm  thick 
barrier  liner,  it  is  shown  that  about  half  of  the  potentially  impinging 
corrosive  species  can  be  kept  away  from  the  metal  substrate  when  the 
diffusion  coefficient  of  the  liner  is  one-tenth  that  of  the  decomposing 
matrix. 

In  Section  IV,  experimental  results  for  the  diffusion  of  gaseous  HF  in 
a  filled  or  unfilled  polyurethane  liner  are  reported.  It  is  shown  that 
SrCrCfy  or  carbon  black  can  effect  a  tenfold  or  greater  reduction  in  the 
polyurethane  diffusion  coefficient.  Thus,  it  seems  reasonable  to  expect 
filled  polymeric  liners  to  prevent  about  50%  of  the  potentially  impinging, 
corrosive  particles  from  reaching  the  metal  substrate. 
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II.  INITIAL  RESULTS  FOR  GASEOUS  HF  DIFFUSION 


A.  BACKGROUND 

This  investigation  considers  the  case  in  which  adhesion  must  be 
maintained  over  an  approximately  10-year  lifetime,  during  which  there  is  a 
slow  evolution  of  corrosive  HF  gas^  from  a  material  layer  that  has,  on  one 
side,  a  sensitive  adhesive /metal  interface  and,  on  the  other  side,  a  layer 
of  material  (or  a  free  surface)  that  is  not  sensitive.  The  polymer 
adhesive  is  chosen  to  be  the  polyurethane,  Uralane  5753.  It  has  been  often 
informally  proposed  that  filling  the  polymer  with  materials  such  as 
alumina,  carbon,  strontium  chromate,  magnesium  silicate,  zinc  oxide,  or 
titanium  dioxide  would  create  a  barrier  liner,  which  would  reduce  the 
amount  of  HF  incident  on  the  wall.  It  is  not  clear  a  priori  if  this  effect 
should  arise  from  chemically  binding  the  HF  within  the  liner  or  from 
kinetically  favoring  diffusion  to  the  material  layer  (or  free  surface)  that 
is  not  sensitive. 

The  Aerospace  Materials  Sciences  Laboratory  (MSL)  has  developed  a 
method  for  quantifying  the  resistance  of  filled  polymers  to  gaseous  HF 
diffusion.  The  method  provides  quantitative  transport  coefficients  for 
filled  polymers  on  the  basis  of  (1)  concentration  profile  measurements  from 
one-dimensional  diffusion  experiments  and  (2)  the  assumption  of  a 
particular  phenomenological  diffusion  model.  The  determination  of  these 
transport  coefficients  is  the  subject  of  this  report.  The  subsequent 
calculation  of  the  HF  flux  incident  on  the  metal  surface  will  be  discussed 
in  the  next  section. 

The  effectiveness  of  filler  particles  needs  to  be  quantified  for  the 
following  reason.  The  matrix  polymer  must  be  selected  for  its  intrinsic 
bondability  to  titanium  alloy  and  its  intrinsic  resistance  to  HF  diffu¬ 
sion.  However,  if  filler  particle  effectiveness  is  well  characterized, 
then  the  matrix  polymer  can  be  selected  with  an  increased  emphasis  on 
intrinsic  bondability  and  with  the  expectation  that  the  filler  will 
adequately  impede  HF  diffusion.  Of  course  this  approach  must  be  used 
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cautiously,  because  it  is  highly  unlikely  that  the  properties  of  highly 
filled  systems  will  be  simple  averages  of  the  intrinsic  properties  of  the 
filler  and  of  the  matrix. 

B.  INITIAL  OBSERVATIONS 

The  search  for  a  suitable  barrier  liner  begins  with  the  selection  of 
an  appropriate  matrix  polymer.  We  initially  experimented  with  polyvinyl 
butyral  (PVB),  which  is  a  well-known,  chemically  resistant  resin  manufac¬ 
tured  by  Monsanto  under  the  name  of  Butvar.  By  varying  the  polymerization 
conditions,  Monsanto  produces  different  kinds  of  Butvar  with  differing 
percentages  of  acetal,  alcohol,  and  acetate  side  groups.  Ue  experimented 
with  the  Monsanto  Metal  Coating  2009,  containing  5.1  wt?  phenolic,  13  wt% 
epoxy,  2  wt£  Butvar  B-90,  and  79. 9t  of  various  solvents. 

In  order  to  judge  the  HP  resistance  of  the  unfilled  Coating  2009,  we 
first  prepared  the  surfaces  of  T1-6A1-4V  coupons  with  sodium  dichromate 
(Na2Cr20y)  and  sulphuric  acid  etching,  followed  by  rinsing  and  vacuum 
drying  (see  Fig.  1-1).  Metal  Coating  2009  was  then  sprayed  on  these 
coupons  to  a  thickness  of  5  ±  2  mils  (see  Fig.  II-1).  The  original  intent 
was  to  use  the  metal  surface  beneath  the  coating  as  a  "detector"  of  the 
time-dependent  corrosion  caused  by  HF  penetration.  This  experiment  was 
actually  performed  by  holding  a  coated  coupon  for  1  day  in  a  Teflon 
accelerated-aging  chamber  of  "pure"  HF  vapor.  (Some  ambient  atmosphere  was 
present  before  the  chamber  was  purged  with  HF.)  In  Fig.  II-2a,  it  is  seen, 
first,  that  debonding  occurred.  Therefore,  the  duration  of  the  originally 
intended  experiment  could  be  no  more  than  a  few  hours.  Second,  a 
nonuniform  buildup  of  corrosion  products  is  seen  in  Fig.  II-2a.  Even  if 
the  film  had  adhered  to  the  coupons  for  a  few  hours,  several  cross  sections 
of  one  coupon  specimen  at  one  aging  time  could  give  widely  varying  results 
for  the  corrosion  buildup.  These  variations  could  be  caused  by  thickness 
variations  of  the  sprayed  films. 
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HF  EXPOSURE 


slrT 


as  Originally  Spray  Coated  with  PVB. 
ight  is  a  sample  holder. 


Fig.  I 1-2 .  (Top  photo,  next  page)  Enlarged  View  of  One  Corner  of  the 

Coated  Coupon  in  Fig.  1 1  —  1  after  1  Day  of  Exposure  to  "Pure" 
(see  text)  HF  Vapor. 

(Bottom  photo,  next  page)  Alumina-Filled  (top)  and  Unfilled 
(bottom)  Uralane  5753  Samples  Cut  from  Larger  Slabs  of 
Material  and  Subjected  to  One-Dimensional  Diffusion  of  HF 
Vapor  (from  the  left).  The  top  sample  is  approximately  1/8 
in.  thick. 
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Third,  and  most  important,  was  that  as  corrosion  products  accumulate 
at  the  metal  interface  within  a  coated  coupon,  the  boundary  conditions  for 
diffusion  change.  Diffusion  is  governed  by  a  partial  differential  equation 
whose  solutions  are  codetermined  by  the  boundary  conditions  at  the  metal 
interface.  In  this  initial  experimental  arrangement,  the  initial  boundary 
condition  at  the  wall  caused  every  impinging  HF  particle  to  interact  with 
the  metal  and  thus  disappear  from  the  diffusion  process.  After  a  suffi¬ 
cient  buildup  of  corrosion  products,  as  shown  in  Fig.  1-2,  the  boundary 
condition  resulted  in  no  impinging  HF  reaction.  In  the  initial  experiment, 
it  would  have  been  difficult  to  separate  this  effect  from  the  intrinsic 
property  of  a  barrier  liner  material  to  impede  diffusion. 

Because  of  these  problems,  we  improved  the  experimental  arrangement. 
Our  goal  was  to  age  the  coupons  in  pure  HF  to  accelerate  diffusion  and  to 
obtain  experimentally  manageable  measurement  times.  To  realize  this  goal, 
we  had  to  find  a  measurement  configuration  without  time-dependent  boundary 
conditions,  so  that  the  transport  properties  of  the  liner  material  could  be 
clearly  seen. 

The  improved  arrangement  selected  for  these  diffusion  experiments 
included  a  large  slug  of  filled  or  unfilled  material  cured  in  a  poly¬ 
ethylene  cup  of  about  1  in.  diam  and  1  in.  depth.  The  polyethylene,  being 
relatively  impervious  (because  of  its  crystalline  content)  and  unreactive 
to  HF,  limited  HF  diffusion  to  one  circular  face  of  the  sample  when  the  cup 
was  placed  in  a  Teflon  Jar  containing  HF.  After  the  slug  was  exposed  to 
HF,  it  was  sectioned  so  that  ion  mass  microprobe  analysis  (IMMA)  or  X-ray 
photoelectron  spectroscopy  (XPS)  could  be  used  to  determine  fluorine 
content  (normalized  by  carbon  content)  vs  position.  This  experimental 
arrangement  was  preferable  because  in  any  of  our  possible  experimental 
times,  diffusion  of  HF  to  the  back  wall  of  the  sample  container  was 
impossible.  Thus,  the  assumed  boundary  condition  at  the  back  wall  was  much 
less  important  when  evaluating  the  concentration  profile. 
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It  would  have  been  advantageous  to  test  filler/resin  composites  where 
the  resin  is  the  liner  material.  However,  this  is  not  possible,  because 
the  sprayable  liner  material  contains  a  high  percentage  (8QJ  in  the  case  of 
Monsanto  2009)  of  solvent.  Attempts  to  cure  such  material  in  bulk  always 
yield  large  air  bubbles.  Attempts  to  cure  such  material  layer  by  layer  are 
time  consuming  and  result  in  stratified  samples  with  unknown  interfacial 
resistance  to  diffusion  between  strata.  Therefore,  a  model  polymer  was 
used  as  the  matrix,  with  and  without  several  filler  types.  The  purpose  of 
the  experiment  is  to  discover  (1)  how  the  matrix-polymer  diffusion 
coefficient  is  affected  (compared  to  the  unfilled  matrix)  by  the  presence 
of  the  filler  particles  and  (2)  to  what  extent  the  filler  particles 
additionally  contribute  to  the  composite's  diffusion  coefficient  (e.g.,  by 
surface  reactions  on  the  filler  particles  that  chemically  trap  some  of  the 
fluorine).  When  this  knowledge  is  subsequently  applied  to  real  liner 
systems,  the  range^  of  values  for  the  unfilled  matrix  polymer  can  be 

Q  p 

assumed  to  be  10“°  to  10“°  cnr/sec,  according  to  recorded  experience  (see 
Table  11-1)  with  small  gas  molecules  in  polymer  matrices.  Using  this  range 
and  experimental  information,  we  can  calculate  the  HF  flux  on  the  metal 
substrate  and  thereby  estimate  the  potential  for  corrosion. 

Figure  II-2b  exhibits  specimens  prepared  from  alumina-filled  and 
unfilled  slugs  of  Uralane  5753  (polyurethane)  and  exposed  to  HF  vapor  for 
12  days.  These  specimens  were  made  several  years  previously  by  MSL  in 
conformance  with  a  specification  requiring  a  10  ±  2  ym  particle  size  (50  um 
maximum).  The  alumina  was  fired  at  700°C  for  10  min  and  stored,  covered, 
at  100-130°C  for  24  h  before  use.  Figure  I I -2b  indicates  that  the  filled 
system  does  retard  HF  diffusion  (shown  by  the  darkest  area  of  HF  reaction). 
An  experimental  technique  was  used  to  quantify  the  different  HF  concen¬ 
tration  profiles  in  the  two  samples  shown  in  Fig.  II-2b.  However,  before 
describing  this  technique,  we  will  construct  a  physical  model  for  guidance 
in  the  experimental  design  and  in  the  interpretation  of  the  resulting 
profiles . 
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Table  II— 1 .  Diffusion  Coefficients®  for  Gas  Holecules  in 
Solid  Polymers  at  25°Cb 


Polymer 

h2 

(MW=2)e 

He 

(MW=4)C 

Gas 

h2o 

(MW=18)C 

co2 

(MW=44)C 

poly ( 1 , 3-butad iene ) 

9.6 

1.05 

poly(butadlene- 

co-acry Ion i tr i 1 e ) 

6.43 

15.5 

— 

0.425 

poly(isoprene) 

— 

21.6 

— 

1.25 

poly (ethylene) 

0=0.914 

— 

6.8 

— 

0.37 

poly (ethylene) 

0=0.964 

— 

3.07 

— 

0.124 

poly (ethyl  methacrylate) 

— 

42.3 

0.0989 

0.033 

poly (styrene) 

— 

10.4 

0.14 

0.058 

poly(tetrafluoro 

ethylene) 

— 

— 

— 

0.095 

Teflon  FEP 

-- 

— 

— 

0. 105 

poly(vinyl  chloride) 

2.8 

0.0238 

0.0025 

®Units  are  10“^  cm2/sec. 
bSee  Ref.  5  <p.  II 1-232) . 
cMolecular  weight  in  grams  per  mole. 
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PHYSICAL  MODEL  FOR  REACTIVE  GAS  DIFFUSION  IN  FILLED  POLYMERS 


Consider  a  one-dimensional  slab  of  material  between  positions  x  =  0 
and  x  =  ft.  The  slab  initially  contains  zero  HF.  The  concentration  of  HF 
(number  of  HF  particles  per  length  of  material)  within  the  slab  is  denoted 
by  e(x,t),  where  t  is  time.  In  this  continuum  model,  any  filler  particles 
(micron  sized)  are  considered  to  be  much  smaller  than  the  slab 
thickness,  ft.  (The  slab  thickness  is  mil  sized,  i.e.,  about  25  times 
larger  than  the  particle  size).  A  constant  concentration,  c0,  of  gas  is 
present  at  the  slab  surface,  x  =  0,  and  no  gas  transport  is  assumed  to 
occur  at  x  =  ft.  However,  as  shown  in  the  preceding  section,  this  second 
boundary  condition  is  not  critical  for  short  experiment  times.  It  is 
required  to  determine  c(x,t)  in  order  to  calculate  the  number  of  HF 
particles  that  impinge  on  the  back  wall  of  the  sample. 

The  analysis  of  unreactive  gas  diffusion  proceeds  in  terms  of  a 
diffusion  coefficient,  D,  for  the  system.  The  analysis  must  be  refined  for 
the  present  case  of  reactive  gas.  We  hypothesize  that  the  following 
boundary -value  equations  are  applicable  to  this  problem,  where  H  is  the 
reactivity  of  the  gas  with  the  system  and  -  He  is  the  rate  of  HF  particle 
reduction  resulting  from  this  reaction: 


Da^c  -  He  =  afcc 


for  0  <  x  <  ft,  t  >  0  and  D  and  H  constant 


Initial  condition: 


Boundary  conditions: 


c(x,o)  =  0 


c(o,t)  =  c 


3  c(  ft ,  t )  =  0 


In  this  boundary  value  problem,  the  units  of  D  are  evidently  area  per  time, 
and  the  units  of  H  are  inverse  time.  The  feature  that  allows  for  the 
reactive  nature  of  the  gas  is  the  term  -  He. 

There  are  at  least  three  possible  physical  interpretations  for  the 
term  -  He.  [Of  course,  the  phenomenological  suitability  of  the  model  can 
only  be  established  by  matching  calculated  and  measured  profiles. 

Acceptance  or  rejection  of  the  microscopic  mechanisms(s)  underlying  this 
problem  must  be  determined  by  other  methods.]  First,  -  He  could  be  a  sink 
term  for  the  annihilation  of  diffusing  molecules  at  a  rate  proportional  to 
the  concentration  of  HF.  A  sink  term  is  conmon  in  heat  conduction  problems 
but  would  conflict  with  our  ideas  about  the  conservation  of  matter. 

Second,  the  -  He  term  could  refer  to  a  rate  of  trapping  (chemically 
inmobilizing)  the  diffusing  particles  by  reactive  sites  in  the  system. 

This  idea  seems  attractive  but  does  not  allow  for  the  saturation  of 
available  trap  sites  at  very  long  times.  Third,  the  -  He  term  could  refer 
to  the  rate  of  diffusion-channel  clogging  by  the  ongoing  process  of 
chemical  reaction  and  physical  adsorption  of  the  diffusing  particles.  This 
interpretation  can  be  considered  mathematically  to  be  caused  by  a  reverse 
current  of  rejected  particles.  This  current  reduces  the  net  flux  of 
particles  into  regions  of  lower  concentration.  Alternately,  and  most 
generally,  this  third  case  can  be  interpreted  as  due  to  a  time-dependent 
diffusion  coefficient  resulting  from  all  physical  and  chemical  processes  of 
diffusant  immobilization.  Consider  the  case  of  no  reaction: 

-V-J  =  3tc 
t 

c(x,t)  =  -  J  dt'  v-J(x,t') 
o 

where  J  is  the  particle  flux.  When  reaction  occurs 


-’■J,rr  =  -’•J  '  Hc  =  V 

(6a) 

t 

c( x , t )  =  -  J  dt'  7-J  ff 

0 

(6b) 

(5a) 

(5b) 
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where  Jeff  is  the  net  flux  into  regions  of  lower  concentration.  Formally 
combining  Eqs.  (6a)  and  (6b),  we  obtain 


t 

-V-J  ♦  H  J  dt'  V*Jeff  =  3tc  (7a) 

o 

Assuming  interchangeability  of  integration  and  differentiation  for  non- 
pathological  functions 

t 

-V-(J  -  H  J  dt'  Jeff)  =  8tc  (7b) 

o 

Comparing  Eqs.  (6a)  and  (7b) 

t 

Jeff(x,tO  =  J(x,t)  -  H  J  dt'  Jeff(x,t')  (7c) 

o 

We  define  an  effective  diffusion  coefficient,  De^f(x,t),  by 

Jeff(X’t)  =  -Deff(x,t)v  c(x,t)  (8a) 

Therefore,  whenever  vc  is  nonzero,  Eqs.  (7c)  and  (8a)  yield 


DeffU,t>  =  D  -  H  j‘  dt-  Deff(x,t') 


(8b) 


I 
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This  equation  is  of  the  form 


Deff(t)  =  D  -  H  /  dt'  Deff(t') 
o 


(8c 


which  is  a  Vol terra  integral  equation  with  the  kernel  K(t',t).^  Dg^j.(t)  i; 
a  function  whose  value  clearly  decreases  monotonically  from  the  value  D  at 
t  =  0.  Before  Eq.  (8c)  can  be  solved,  c(x,t)  must  be  found. 

The  solution  c(x,t)  to  the  single  slab  problem  specified  by 
Eqs.  (1)-(k)  could  be  obtained  from  a  Sturm-Liouville  eigenvalue 
approach.  However,  the  calculation  of  HF  flux  at  the  metal  substrate  (to 
be  presented  in  Section  III)  clearly  must  allow  for  two  adjacent  layers 
(matrix  and  liner).  Churchill^  has  shown  that  the  solution  to  such  a 
problem,  in  which  there  is  more  than  one  layer,  cannot  be  constructed  with 
the  Sturm-Liouville  approach.  Therefore,  the  alternate  method  of  Ref.  8  is 
used  here  as  a  deliberate  choice  and  in  subsequent  multislab  calculations 
as  a  necessity. 

Laplace-transforming^  Eqs.  (1)  and  (2),  we  obtain 


2-  2- 
3  c  -  w  c  =  0 
x 


2  H  +  s 
u  =  — =: — 


(9a 

(9b 


where  c(s)  is  the  Laplace  transform  of  c(t),  suppressing  the  x-vanable 
notation.  Consequently 


c  =  A  cosh  wx  B  sinh  u>x  (9c 

where  A  and  B  are  constants.  Applying  Eqs.  (3)  and  (4),  we  obtain 


c(s) 


Q 

o  fCosh[b)i.(1  -  x/i)]i 
s  ‘  cosh  ‘ 


Letting  the  denominator  of  c(s)  be  denoted  by  f(s),  we  obtain 


f(s)  =  s  cosh 


(10a) 


Note  that  f(s)  >0  for  all  s  >  0.  For 
|s|  <  H,  f(s)  <  0.  For  s  <  0  and  |s| 

s  =  0,  f(s)  =  0.  For  s  <  0  and 
>  H,  f(s)  =  0  whenever 

_ 

0  =  sn  cosh[i^ 

/-(H  ♦  sn) 

(10b) 

'  D  M 

or 

-- 

sn  *  -H  - 

-  D* 

T 

(10c) 

n 

=  1,2,3,... 

T 

Dx2 

Therefore,  the  pole  structure  (which  is  all  first  order)  of  c(s)  is 

lm(s)  jjs 

,  .  .  V  .  .  —  V-  y  y  < 

•  •  •  A  A  A  i 

Sn  s3  3 2  31 

Re(s) 

Expanding  c(s)  according  to  the  Mittag-Leffler  theorem,10  we  have 

— 

c(s)  -  Res[c(s),  Q] 

®  Res[c(s),s  ] 

-  +  l  - — 

L.  s  -  s 

n=1  n 

Ola) 
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where  Res[c(s),  s']  denotes  the  residue  of  c(s)  at  s  =  s' 
residues,  we  have 


Evaluating  the 


Res[c(s),  o]  = 


c  cosh[u  t(1  -  x/t)] 


cosh  oi  t 
o 


u  =  /H/D 
o 


Res[c(s),sn] 


Ms~Th) 

-t- 


(H*sn>  , 
— n -  X1 


Laplace- inverting  Eq.  (11a),  we  obtain 

cosh[w_t(1  -  x/l)]  ^ ®  /-(s_  ♦  H) 


/-(s_  ♦  H) 


cosh  u  l 

i  o 

1  n=1 

s  -‘"i/  1 

n 

wo 

=  /H/D 

-  u 

(2n  -  I)2 

n  =  1,2,3, .. . 

sn 

-  —  fl  — 

T 

t 

Dt2 

COMPUTER  CALCULATIONS 

AND  THEIR  PHYSICAL  INTERPRETATION 

The  results  shown  in  Eq.  (12)  have  been  programned  on  the  Aerospace 
Space  Science  Laboratory's  (SSL's)  "DIRAC  VAX.  The  normalized  profiles, 
(c(x,t)/c0),  are  given  in  Figs.  1 1 -3  through  1 1 -7 .  Figure  1 1 -3  represents 
the  case  H  =  0,  which  is  diffusion  without  chemical  reaction,  diffusion- 
channel  clogging,  and  so  on.  In  this  report,  we  use  megaseconds  as  the 


C(X)/C0  C(X)/C 


1.02 


2  Msec 


X  (cm) 

Fig.  II-4.  Calculated  Diffusion  Profiles  for  H  =  1  Msec-1 
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Fig.  1 1-5.  Calculated  Diffusion  Profiles  for  H  =  10  Msec 


unit  of  time  for  the  calculations.  One  megasecond  is  278  h.  At  long 
times,  t  >  2  Msec  >  t  =  0.406  Msec,  the  profile  is  flat,  and  diffusion  is 
complete. 

Numerical  values  of  the  parameters  used  in  the  calculation  of  the 
normalized  profiles  are  as  follows: 

D  =  1  cm2 /Msec 

i  s  1  cm 

T  =  Hi2 /Du2  =  0.406  Msec 
t  in  Msec  (1  Msec  =  278  h) 
x  in  cm 

Profiles  were  calculated  at  the  following  times:  10'^,  5  *  10'^,  10'2, 

5  x  10'2,  10'1,  0.5,  1,  and  2  Msec. 

Figure  I I -3  shows  calculated  profiles  for  H  =  0  Msec-1;  Fig.  II-4,  for 
H  =  1  Msec'1;  Fig.  I 1-5,  for  H  =  10  Msec'1;  Fig.  11-6,  for  H  =  100  Msec'1, 
and  Fig.  11-7  for  H  =  1000  Msec'1.  In  Fig.  II-4  (where  H  =  1),  we  see  that 
the  0.5,  1,  and  2  Msec  results  are  significantly  shifted.  In  fact,  even  as 
t  ♦  ®,  the  flat  profile  is  not  attained.  This  is  shown  by  taking  the 
derivative,  oQ,  of  Eq.  (12)  at  x  =  0  and  letting  t  tend  to  infinity: 


o 


o 


lim  —  3c 
.  c  x  1 
t-»=>  o  x=o 


=  -  /^~  tanh 


(13a) 


Because  oQ  <  0  for  H  >  0,  the  flat  profile  is  never  attained.  The  final 
value,  PQ,  for  the  plateau  level  is 

PQ  =  lim  c(l,t)  =  sech  {/ ^  l)  (13b) 

t-»®  o 


Some  numerical  results  for  o  and  P  are  given  in  Table  I I -2. 

o  o 
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Table  II-2.  Initial  Slopes  (a  )  and  Final  Plateau 
Levels  (P  )  for  tRe  Normalized  Profiles 
for  Several  Values  of  the  Interaction 
Parameter  H 


H 

(Msec"1 ) 

"o, 

(cm"1) 

P 

0 

(dimensionless) 

0 

0 

1 

1 

-0.762 

0.65 

10 

-3.15 

0.085 

100 

-10.0 

o 

1 

Jr 

1000 

-31.6 

.= r 

1 

o 

This  suppression  of  the  infinite-time  plateau,  shown  in  Figs.  I 1-4 
through  II -7 ,  is  consistent  with  a  "rejection  current"  resulting  from 
annihilation,  trapping,  clogging,  and  so  on: 

l 

L  =  -  J  dx  H  c(x,t)  ( 14a) 

o 

G  =  J(o, t)  =  -  DV  c( x , t )  |  (14b) 

x=o 


where  L  is  the  rate  of  reduction  of  mobile  HF  particles  throughout  the  slab 
and  G  is  the  rate  of  increase  of  HF  particles  caused  by  diffusion  from  the 
environment  into  the  slab  at  position  x  =  o.  As  the  time  tends  to 
infinity,  the  integral  and  the  derivative  of  the  infinite  series  terms 
{compare  Eqs.  (12),  (14a),  and  (14b)]  are  zero,  so  that  L  =  -  /HD  cQ  and 
G  =  /HD  cq  (particles/Msec) . 
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In  Fig.  I 1-5  (H  =  10),  we  see  that  the  results  for  t  =  0.5,  1,  and 
2  Msec  actually  superimpose.  In  Fig.  II-6  (H  =  100),  the  last  five  times 
have  coincident  profiles.  In  Fig.  I 1-7  (H  =  1000),  this  saturation  is 
reached  at  t  =  5  *  10"^  Msec,  which  is  only  about  1.4  h.  Clearly,  HF 
diffusion  can  be  significantly  restricted  for  a  sufficiently  large 
interaction  parameter,  H.  If  measurements  are  made  at  a  time  close  to  the 
system  relaxation  time,  t,  and  if  a  result  is  obtained  like  that  for  t  = 

0.5  Msec  (1.25  t)  in  Fig.  I I -3,  then  H  is  undoubtedly  very  small.  If  the 
results  are  like  those  from  t  =  0.005  to  2  Msec  in  Fig.  11-7,  then  H  is 
very  large.  We  need  a  method  for  extracting  D  and  H  from  measured 
profiles . 

We  know^  that  one-dimensional  diffusion  into  an  unbounded  medium  is 
governed  by  <x>  =  2  /Dt,  where  <x>  is  an  "effective"  penetration  depth.  In 
Fig.  1 1 — 3 ,  let  <x>  be  the  position  at  which  c(x,t)/cQ  assumes  the  value 
0.17.  The  values  in  Fig.  I I -3  of  <x>  for  t  =  10”^  to  10_1  Msec  are 
unaffected  by  the  presence  of  the  back  wall  of  the  sample  container.  These 
values  are  listed  in  Table  I 1-3  and  Fig.  11-8.  On  the  other  hand,  the  role 

Table  II-3-  Propagation  Into  an  Unbounded  Medium3 


<  X  > 

(cm) 

t 

(Msec) 

✓  t 

(Msec1 /2) 

<x>/(2/t) 

(cm/Msec1^) 

5 

p 

(cm  /Msec) 

0.062 

0.001 

0.0316 

0.98 

0.96 

0.14 

0.005 

0.0705 

0.99 

0.98 

0.20 

0.010 

0.1000 

1 

1 

0.44 

0.050 

0.2236 

0.98 

0.96 

0.6 1 

0.100 

0.3162 

0.96 

0.92 

aData  from  Fig.  I 1-3  (H  =  0)  for  t  <  0. 1  Msec: 

<x>  defined  at  c(x,t)/cQ  =  0.17 
D  =  1  cm^/Msec 

D  =  apparent  (secant)  value  of  D  (see  Fig.  I I -8 ) 


I 1—8 .  "Effective”  Diffusion  Front  (H  =  0)  Advances  According  to  <x> 
2/Dt  into  an  Unbounded  Medium  (for  t  <  t,  i.e.,  when  the  back 
wall  plays  no  role).  The  plotted  points  are  obtained  from 
Fig.  1 1-3  (see  text) . 


of  H  =  10  in  Fig.  I I —5  is  seen  in  Table  II-4  and  Fig.  I I -9  as  a  deviation 
from  the  "classical"  <x>  =  2  /Dt  trajectory.  The  deviation  of  <x>  from 
2  /Dt  is  most  pronounced  in  the  case  of  H  =  1000  Msec-1  (see  Table  I 1-5  and 
Fig.  11-10). 

Table  II-4.  Propagation  into  an  Unbounded  Medium* 


<x> 

(cm) 

t 

(Msec) 

✓t 

(Msec1/2) 

<x>/(2/t) 
(cm/Msec 1/2 ) 

B 

(cm2/Msec) 

0.062 

0.001 

0.0316 

0.98 

0.96 

0.14 

0.005 

0.0705 

0.99 

0.98 

0.19 

0.010 

0.1000 

0.95 

0.90 

0.39 

0.050 

0.2236 

0.87 

0.76 

0.49 

0.100 

0.3162 

0.77 

0.59 

0.585 

0.5 

0.705 

0.41 

0.17 

aData  from  Fig.  I 1-5  (H  =  0)  for  t  <  0.5  Msec: 

<x>  defined  at  c(x,t)/cQ  =  0.17 
D  =  1.0  cm2 /Msec 

B  =  apparent  (secant)  value  of  D  (see  Fig.  I I -9 ) 

Compare  Figs.  11-8,  1 1-9,  and  11-10,  where  <x>  is  plotted  vs  /t.  In 
Fig.  I 1—8 ,  when  several  measurements  lie  on  a  straight  line  through  the 
origin,  we  may  identify  the  slope  with  2  /D.  From  Fig.  11-8,  by  itself,  we 
can  only  conclude  that  the  measurement  times  are  too  short  for  any  H  value 
to  be  determined  (although  we  know  that  H  =  0  in  this  particular  case). 
Figure  11-10  is  at  the  other  extreme,  where  there  are  not  at  least  two 
nonzero  <x>  measurements  on  a  straight  line  through  the  origin.  We  can 
only  conclude  that  the  measurement  times  were  too  long  to  determine  D. 
Figure  I 1-9  is  the  intermediate  case,  where  D  can  be  determined  from  the 
slope  of  the  line  that  passes  through  the  first  two  nonzero  <x>  points  as 
well  as  the  origin.  The  value  of  H  in  Fig.  I I -9  can  be  determined  as 
follows. 
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Fig. 


Fig.  I 
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0  0  01  0  04  0  09  0  16 
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I 1-9.  "Effective”  Diffusion  Front  (H  =  10)  Advances  According  to 
<x>  =  2V Dt  Until  Saturation_(due  to  the  nonzero  H)  Sets  In. 
Here,  deviation  starts  at  /t  =  0.223  Msec1'2.  The  plotted 
points  are  obtained  from  Fig.  I I -5  (see  text). 


t  (Msec) 


1-10.  "Effective"  Diffusion  Front  of  a  Strongly  Interacting  System 
(H  =  1000)  Does  Not  Exhibit  Classical  <x>  =  2/Dt  Behavior  in 
This  Region.  For  times  sufficiently  smaller  than  10“  ^  Msec, 
classical  behavior  should  be  regained.  The  plotted  points  are 
obtained  from  Fig.  1 1 -7  (see  text). 
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Table  I I -5 .  Propagation  Into  an  Unbounded  Medium3 


<x> 

(cm) 

t 

(Msec) 

/t 

(Msec172) 

<x>/(2/t) 

(cm/Msec172) 

6 

(cm2/Msec) 

0.047 

0.001 

0.0316 

0.74 

0.41 

0.057 

0.005 

0.0705 

0.40 

0.16 

0.057 

0.010 

0.1000 

0.28 

0.078 

0.057 

0.050 

0.2230 

0.13 

0.017 

0.057 

0.100 

0.3162 

0.09 

0.008 

aData 

from  Fig.  II- 

-7  (H  =  0)  for  t 

<  0. 1  Msec: 

<x> 

■  defined  at  c(x,t)/cQ  = 

0.17 

D 

=  1.0  cm2/Msec 

0 

=  apparent  (secant)  value  of 

D  (see  Fig.  II 

-10) 

Ordinary  (H  =  0)  diffusion  into  a  slab  may  be  considered  as  a  relax- 

2  2 

|  ation  of  concentration  with  a  time  constant  t  =  4l  /Du  .  This  relaxation 

is  typical  of  a  process  occurring  at  a  rate  v,  where 

v  =  &  (15a) 

/t 

I 

t  _ 

< x >  =  f  v{t')dt'  =  2  /Dt  (15b) 

o 

^  On  the  other  hand,  consider  the  data  from  Figs.  1 1 —3  through  1 1 — 7 ,  listed 

in  Table  1 1 -6 .  Apparently,  a  second  relaxation  is  occurring  with  a  time 
constant  H*1.  We  hypothesize,  therefore,  that 


I 

v 


! 


=  ~  exp  ( -Ht) 


(16a) 
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t 

<x>  =  J  v(t')dt' 
o 


erf  (/Ht) 


(16b) 


Using  the  expansion^  of  the  error  function  as  well  as  its  asymptotic 
value,  we  find 


2  /Dt  [l  - 

Ht  ^  jHtj!  _  , 

3  10  •  J 

(16c) 

< x  >  =  2  /Dt 

Ht  <<  1 

( l6d) 

✓  s  /sir 

<x>  =  /r 

t  ♦  » 

(I6e) 

Letting  f  =  <x>/(2  /Dt)  and  using  Eq.  (16c)  (up  to  the  quadratic  term),  we 
obtain  the  value  of  H  from  the  short-time  deviation  of  <x>  from  2  /Dt  as 
follows: 


H 


0.4  (1  -  f) 
0~2 


(17) 


Table  1 1-6.  Correlation  of  Profile  Saturation  with  the 
Value  of  the  Interaction  Parameter,  Ha 


H  . 

(Msec  ) 

Approximate  saturation  time 
(Msec) 

1/H 

(Msec) 

0 

OD 

OD 

1 

2 

1 

10 

0.5 

0.1 

100 

0.05 

0.01 

1000 

0.005 

0.001 

aData  from  Figs.  1 1 - 3  through  1 1 -7  - 
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Applying  Eqs.  (16)-(17)  in  Fig.  11-9,  we  find  that  classical  <x>  = 

2  /Dt  diffusion  occurs  for  t  <<  H’1  =  0.1  Msec.  We  also  find  that 

/Dw/H  =  0.56  cm  is  a  reasonable  extrapolation  for  long-time  behavior.  To 

find  H,  we  interpolate  a  synthetic  data  point  between  t  =  0.01  and  0.05, 

balancing  the  requirements  of  t  <<  0.1  Msec  and  of  <x>  being  measurably 

different  from  2/Dt.  We  use  t  =  0.0256  Msec  and  <x>  =  0.29  cm.  Applying 

Eq.  (17)  with  f  =  0.29/0.32  =  0.906,  we  find  H  =  12  Msec’1,  which  deviates 

significantly  from  the  true  value,  H  =  10  Msec’1.  We  conclude  that  the 

small  deviations  from  the  classical  behavior  will  not  accurately  determine 

H.  This  is  a  strong  conclusion,  because  it  implies  that  diffusion-profile 

measurements  must  be  made  at  times  close  to  the  system  relaxation  time 
2  2 

(-t  =  4l  /Du  )  in  order  to  determine  H.  These  times  should  be  as  follows: 
for  l  =  2.54  cm  and  D  =  1.0  cm^/Msee,  i  =  2.6  Msec  =  727  h  =  30.3  days; 
for  1  =  2.54  cm  and  D  =  0.01  cm^/Msec,  t  =  3030  days  =  8.3  years.  However, 
such  measurement  times  are  impracticable.  Therefore,  H  values,  such  as 
those  in  Fig.  I I -9,  are  best  estimated  only  by  order  of  magnitude. 

To  accurately  determine  H,  we  must  have  plateau  values  of  <x>  vs  /t. 
For  example,  Fig.  11-11  shows  one  more  data  point  (at  0.5  Msec,  taken  from 
Fig.  1 1-4 )  than  Fig.  1 1-9 .  The  plateau  value  is  0.585  cm.  Using 
Eq.  (I6e),  H  =  9.2  Msec’1.  This  determination  of  H  is  not  in  agreement 
with  the  true  value  of  10;  however,  it  is  at  least  an  improvement  over  the 
estimate  from  short-time  behavior.  The  presence  of  the  back  wall  of  the 
sample  may  be  affecting  this  result.  The  agreement  could  possibly  be 
improved  if  a  refinement  of  the  simple  postulate,  Eq.  (16a),  could  be 
obtained.  Considering  Table  11-5  and  Fig.  11-10,  we  note  that  the 

2  i 

saturation  <x>  =  0.057  cm  is  consistent  with  H  =  D*/<x>  =  967  Msec  , 

which  is  close  to  the  true  value,  H  =  1000  Msec"1.  For  such  a  strong 
interaction  parameter,  the  measurements  may  be  done  in  a  time  much  less 
than  t.  The  problem,  as  noted  in  the  discussion  of  Fig.  11-10,  is  that  D 
cannot  be  determined  from  this  data.  Shorter  measurement  times  would  have 
to  be  used  to  obtain  the  slope  of  <x>  vs  /t  in  the  classical  regime. 
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<x>  (cm) 


In  summary,  to  determine  the  D  and  H  values  for  a  system,  data  must  be 
obtained  in  the  form  shown  in  Fig.  11-11.  In  this  figure,  the  initial 
slope  provides  D,  and  the  final  plateau  provides  /Di/H.  If  we  cannot 
obtain  long  enough  measurement  times  to  reach  the  plateau,  we  must  either 
assume  H  =  0  or  at  best  estimate  the  order  of  magnitude  of  H.  Of  course, 
if  the  results  remain  on  the  <x>  =  2  /Dt  trajectory  (until  coupling  to  the 
back  wall  of  the  sample  invalidates  the  unbounded-medium  condition),  then  H 
=  0,  or  at  least  H  <<  t”\ 

D  and  H  must  be  determined  for  each  polymer/gas  system.  H  specifies 
the  essential  trapping/reactive  nature  of  the  filler  particles  (in  conjunc¬ 
tion  with  the  matrix).  If  the  H  parameter  is  unmeasurably  small,  then 
there  is  no  evidence  to  support  the  surface  reactivity  interpretation  of 
the  filler  particle  effect.  But  considering  the  visual  impression  of 
Fig.  II-2b,  it  is  anticipated  that  filler  particles  make  a  real  differ¬ 
ence.  Even  if  there  were  no  "H  effect,"  there  are  important  effects  of  the 
filler  upon  D.  These  effects  have  been  documented  in  the  diffusion  of 
inert  gases,  as  discussed  in  the  next  two  paragraphs. 

In  the  most  widely  accepted  model  for  gas  transport  through 

Ip  ii: 

polymers  at  temperatures  below  the  glass  transition  temperature,  T  , 

© 

each  gas  molecule  must  Jump  between  fixed  free-volume  sites.  Above  T  ,  the 

5 

polymer  segments  are  more  mobile  and  the  available,  mobile  free-volume 

vacancies  for  gas-molecule  hopping  allow  propagation  (diffusion)  to  occur 

at  a  faster  rate.  The  rate  of  increase  of  the  diffusion  coefficient  with 

temperature  is  relatively  low  for  T  <  Tg,  because  the  polymer  lattice 

expands  relatively  little  and  because  no  additional  free  volume  can  enter 

the  polymer  system.  The  rate  of  increase  of  the  diffusion  coefficient  with 

temperature  is  relatively  high  when  T  >  T_,  because  there  is  progressively 

17  * 

more  free  volume  in  the  system  '  as  the  temperature  rises. 

The  presence  of  filler  reduces  the  mobility  of  the  adsorbed  (or 
chemically  bonded)  polymer  segments.  In  effect,  the  glass  transition 
temperature  is  higher  in  a  thin  layer  of  polymer  adjacent  to  each  filler 
particle,  i.e.,  the  polymer  in  this  layer  is  below  its  T  and  allows  less 

o 


gas-molecule  diffusion.  This  layer  of  restricted  mobility1^  can  be  as 

thick  as  5000  to  8000  A.  Because  the  interparticle  spacing1^  of  a  TiC^- 

f tiled  (19.25  wtf)  polymer  has  been  reported  to  be  1250  A,  the  filler  can 

place  the  entire  matrix  in  a  glassy  state,  which  could  not  be  achieved  in 

an  unfilled  system  at  the  same  temperature.  For  example,  the  T  of 

© 

poly(vinyl  acetate  -  vinyl  chloride)  can  be  raised  from  30  to  51 °C  by 

IP 

appropriate  loading.  The  diffusion  coefficient  of  such  a  glassy  system 
will  be  correspondingly  lower.  Thus  the  filler  could  have  a  significant 
effect  on  the  overall  diffusion  even  if  the  filler  particle  surfaces  were 
inert  to  the  gas.  We  suggest  that  this  effect  also  explains  the  restricted 
HF  penetration  into  our  filled  polymers.  Further  testing  is  in  progress  to 
determine  the  Tg  of  filled  and  unfilled  Uralane. 

It  is  important  to  experimentally  check  the  present  mathematical  model 
in  order  to  ascertain  whether  factors  in  addition  to  concentration 
gradients  are  driving  diffusion.  It  is  known,  in  general,  that  chemical 
potential  gradients  drive  diffusion  and  that  concentration  gradients  are  a 
special  case  of  this  effect.  For  example,  low  levels  of  some  fillers  were 
seen  to  increase  the  permeability  of  some  liners  to  liquid  migrants.1®  It 
was  postulated  that  this  increase  resulted  from  the  chemical  reactivity  of 
the  migrant  with  the  filler. 

E.  EXPERIMENTAL  RESULTS 

One-dimensional  diffusion  experiments  were  conducted  on  alumina-f illed 
(M0  wt%)  and  on  unfilled  samples  of  the  polyurethane  Uralane  5753.  Because 
the  Uralane  was  cured  in  HF-impervious  polyethylene  cups,  the  HF  could 
enter  the  Uralane  only  from  its  one  exposed  face.  The  cup  was  placed  in  an 
HF  aging  chamber,  which  was  a  larger  Teflon  jar  that  had  a  threaded  lid 
sealed  with  an  0-ring.  Inlet  and  outlet  bulkhead  connectors  made  from 
Teflon  were  tightly  installed  in  the  Jar  wall.  Teflon  tubes  connected  the 
inlet  to  an  HF  lecture  bottle  and  the  outlet  to  a  water  bath  that  trapped 
excess  HF.  Teflon  valves  in  the  inlet  and  outlet  lines  maintained  a  tight 
seal  after  HF  was  introduced  into  the  aging  chamber.  An  ambient  laboratory 
atmosphere  (including  humidity)  was  present  in  the  aging  chamber  before  the 


chamber  was  purged  for  5  min  with  HF  to  establish  an  HF  atmosphere  at  room 
temperature.  The  HF  atmosphere  was  "recharged"  once  per  day  during  the 
aging  period.  The  differences  between  HF  diffusion  in  filled  and  unfilled 
samples  are  valid  because  each  sample  was  in  the  identical  environment. 

One  cup  of  unfilled  Uralane  and  one  cup  of  alumina-filled  Uralane  were 
simultaneously  exposed  in  the  same  Teflon  jar  for  21  h.  Slices  were 
subsequently  cut  from  the  inside  of  the  exposed  samples,  with  the  cuts 
perpendicular  to  the  direction  of  diffusion.  These  samples  were  subjected 
to  IMMA  and  were  normalized  by  carbon  concentration.  The  experimental 
diffusion  profiles  (Figs.  II-12a  and  II-12b)  deviated  considerably  from  the 
calculated  diffusion  profiles  (Figs.  11-3  to  I 1-7 ) .  This  deviation  was 
attributed  to  the  surface  roughness  of  the  samples,  which  could  not  be 
polished  because  of  the  "smearability"  of  the  profiles.  However,  the 
experimental  results  for  the  unfilled  sample  were  consistent  with  a 
calculated  profile  of  <x>  =  2  mm  (a  position  where  there  is  about  17%  of 
the  estimated  1.5  initial  concentration  in  Fig.  ll-12a). 

The  IMMA  results  (Fig.  11-1 2b )  for  the  filled  sample  were  consistent 
with  the  superposition  of  two  effects.  First,  there  was  a  sharper  decrease 
in  HF,  of  the  type  seen  in  Figs.  I 1—3  to  II-7,  with  <x>  =  1  mm.  Second, 
there  was  an  additional  constant  plateau  level  of  about  20 %  of  the  initial 
(=  1.5)  HF  concentration,  which  extends  relatively  far  into  the  sample. 

This  is  clearly  an  independent  process  that  has  been  completed  (level 
profile)  but  that  is  not  present  uniformly  throughout  the  sample  (only  20% 
of  the  surface  concentration  was  obtained).  One  plausible  explanation  for 
such  a  process  is  the  HF-induced  debonding  of  polymer  from  filler 
particles.  The  resulting  crevices  could  rapidly  fill  with  HF  gas  but  are 
limited  in  their  volume  fraction.  This  effect  could  be  modeled  by  parallel 
diffusion  pathways,  where  one  diffusion  coefficient  is  at  least  one  order 
of  magnitude  faster  than  that  of  the  gas  molecule  in  the  polymer  matrix. 
Measurements  at  extremely  short  times  would  have  to  be  conducted  to 
ascertain  exactly  how  many  orders  of  magnitude  are  involved.  This  effect 
(of  unestablished  origin)  is  clearly  harmful  for  a  film  proposed  as  a 
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Fig. 


1-12.  (a)  Experimental  (ion  microprobe)  Diffusion  Profile  for 

Unfilled  Uralane  After  21  h  of  Exposure,  (b)  Experimental 
(ion  microprobe)  Diffusion  Profile  for  Alumina-Filled  Uralane 
After  21  h  of  Exposure. 


propellant /motor  case  barrier.  Further  analysis  is  given  in  the  next 
section  in  order  to  determine  exactly  how  harmful  this  effect  is. 

Because  of  the  surface  roughness  problem  observed  in  the  IMMA  data,  it 
was  decided  to  obtain  XPS  data  on  two  similar  samples.  After  12  days  of 
aging,  the  slugs  from  the  similar  samples  were  withdrawn  from  the  aging 
chamber  and  sectioned  for  XPS  analysis.  Initially,  the  oxygen  Is  signal 
was  analyzed,  because  the  HF  degradation  reaction  was  responsible  for 
smearing  this  peak  (intensity  vs  binding  energy),  consistent  with  the 
formation  of  reaction  by-products  with  varied  oxygen  bonding  states.  When 
no  systematic  variation  could  be  found,  carbon-normalized  fluorine  concen¬ 
tration  vs  position  was  analyzed.  The  XPS  results  for  unfilled  and  filled 
Uralane  are  shown  in  Figs.  II-13a  and  II-13b.  The  unfilled  specimen  has 
the  calculated  curve  shape  seen  in  Figs.  I 1—3  to  I 1-7 ,  indicating  that  the 
surface  roughness  problem  is  not  critical  in  XPS.  The  filled  specimen 
(Fig.  1 1  —  1 3b )  shows  the  same  superposition  of  two  effects  seen  in  the  IMMA 
results  (Fig.  II— 12b )  but  with  less  experimental  scatter.  Similarly,  as 
before,  the  estimated  HF  penetration  depths  are  2  mm  (unfilled)  and  1  mm 
(filled).  The  penetration  is  (experimentally)  the  same  even  when  sample 
exposure  increases  from  21  to  288  h.  For  unhindered  diffusion,  the 
diffusion  front  should  have  advanced  by  a  factor  of  3.7  between  21  and 
288  h. 

In  this  discussion,  we  assume  that  the  profiles  in  Figs.  II-12a  and 
II-12b  are  the  same  as  those  in  Figs.  II-13a  and  I I- 13b,  respectively.  (Of 

course,  the  results  shown  in  Fig.  11-12  should  be  repeated  using  XPS.)  If 

this  assumption  is  true,  then  we  have  obtained  the  case  shown  in  Fig.  11-10 
(as  opposed  to  Figs.  I 1-8  and  II-9)  for  both  the  filled  and  unfilled 
samples.  It  follows  that  t  =  21  h  i  H"1,  or  H  *  13  Msec*1,  in  either 
filled  or  unfilled  samples.  It  does  not  follow  that  H  is  the  same  for  both 

filled  and  unfilled  samples,  only  that  the  values  of  H  for  both  samples  are 

consistent  with  H  *  13  Msec"1.  Furthermore,  |oQ|  =  /H/D  =  H  cm"1 
(unfilled)  and  11  cm"1  (filled)  [see  Eq.  (13a)  and  Figs.  II— 13a  and 
I I - 1 3b ] .  The  corresponding  diffusion  coefficients  are  D  *  0.8  cm2/Msec 


(%)  D 


(unfilled)  and  D  i  0.1  cia^/Hsec  (filled).  These  D  values  are  consistent 
with  the  Tg  phenomena  previously  described  and  with  the  range  of  values 
seen  in  Table  11-1.  As  noted,  only  measurements  over  shorter  times  would 
completely  specify  both  D  and  H. 

F.  SUMMARY 

The  Aerospace  Materials  Sciences  Laboratory  has  investigated  the 
assumption  that  the  barrier  effectiveness  of  a  polymeric  liner  will  be 
improved  by  adding  filler  particles  that  can  inhibit  the  flow  of  HF  gas. 

The  transport  coefficients  of  the  liner  materials  must  be  determined  before 
multislab  calculations  can  reveal  the  extent  to  which  the  barrier  reduces 
the  HF  flux  at  the  metal  substrate.  These  transport  coefficients  can  be 
determined  from  one-dimensional  diffusion  experiments  in  which  the  HF 
concentration  profiles  are  measured  by  XPS.  Several  such  profiles  must  be 
measured  in  order  to  determine  both  the  basic  diffusion  coefficient,  D,  and 
a  second  parameter,  H,  specifying  the  degree  of  chemical  reactivity  between 
the  gas  and  the  liner  material.  The  current  experiments  (21  h  and  12  days 
long)  on  a  model  polymer  with  and  without  alumina  particles  (40  wt%  of  10  ± 
2-um  particles)  are  consistent  with  H  *  13  Msec-1  for  either  filled  or 
unfilled  polymer:  D  -  0.8  cm2 /Msec  (unfilled)  and  D  *  0.1  cm2 /Msec 
(filled).  In  Section  IV,  we  will  determine  D  and  H  for  several  different 
filler  types. 

Mathematical  modeling  established  that  diffusion  can  be  significantly 
impeded  by  high  values  of  the  interaction  parameter,  H.  Debonding  between 
filler  and  matrix  may  partially  offset  this  gain.  In  addition,  other 
researchers  have  shown  that  the  presence  of  filler  particles  immobilizes 
polymer  segments,  raises  the  polymeric  glass  transition  temperature,  and 
reduces  the  diffusion  coefficient,  D.  For  example,  the  Tg  of  poly(vinyl 
acetate  -  vinyl  chloride)  can  be  raised  from  30  to  51°C  by  appropriate 
loading.  Any  such  increase  will  lead  to  a  reduced  diffusion  coeffi¬ 
cient.  The  most  dramatic  decrease  occurs  for  a  polymer  system  in  which  the 
loading  elevates  Tg  from  below  the  operating  temperature  to  above  it.  We 
suggest  that  this  effect  explains  the  restricted  HF  penetration  into  our 
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filled  polymers.  Further  testing  is  in  progress  to  determine  the  T„  of 

& 

filled  and  unfilled  Uralane. 

Phenomenologically,  the  description  of  diffusion  in  terms  of  D  and  H 
is  consistent  with  the  existence  of  a  time-dependent  diffusion  coefficient, 
which  may  be  considered  to  arise  microscopically  from  the  clogging  of 
diffusion  channels  by  reaction  products.  The  consequences  of  such  hindered 
diffusion  in  a  layered  system  will  be  described  in  the  next  section. 


III.  CALCULATION  OF  HF  FLUX  AT  THE  METAL  SUBSTRATE 


1 


A.  BACKGROUND 

The  liner  considered  here  must  serve  as  a  barrier  to  corrosive  HF  gas. 
In  this  section,  the  fraction  of  the  evolved  HF  particles  that  can  diffuse 
through  the  liner  and  arrive  at  the  metal  substrate  is  calculated  for  some 
typical  values  of  component  thicknesses,  diffusion  coefficients  (matrix 
and  liner  D2 ) >  and  liner  reactivity  with  HF.  Figure  1-3  may  be  consulted 
in  order  to  visualize  the  geometry  and  parameters.  A  barrier  effective¬ 
ness,  c,  is  defined  such  that  c  =  100^  when  no  HF  ever  reaches  the  metal 
substrate. 

In  Section  I,  a  method  for  quantifying  the  resistance  of  filled  poly¬ 
mers  to  gaseous  HF  diffusion  was  developed  on  the  basis  of  (1)  concentra¬ 
tion  profile  measurements  from  one-dimensional  diffusion  experiments  and 
(2)  the  assumption  of  a  particular  phenomenological  diffusion  model. 
Experiments  on  a  model  polymer  with  and  without  filler  indicated  that  the 
diffusion  coefficient,  D,  is  at  least  0.8  croVMsec  (unfilled)  and  at  least 
0.1  em2/'Msec  (filled).  In  either  case,  the  parameter,  H,  specifying  the 
reactivity  of  the  gas  and  the  decrease  of  D  with  time,  was  found  to  be  at 
least  13  Msec.  Using  these  results,  as  well  as  literature  values  for  the 
diffusion  coefficients  of  common  polymers,  we  calculated  the  barrier 
effectiveness,  t,  for  several  sets  of  parameters.  The  ratio  of  the  total 
number  of  HF  particles  arriving  at  a  lined  metal  substrate  surface  to  the 
total  number  of  HF  particles  arriving  when  unlined  is  (1  -  c/100%). 

B.  PHYSICAL  MODEL  FOR  THE  DIFFUSION  OF  REACTIVE  GASES  OUT  OF 

A  MATRIX.  THROUGH  A  LINER.  AND  ONTO  A  METAL  SUBSTRATE 

Consider  a  one-dimensional  slab  of  a  gas-evolving  matrix  material 
between  positions  x  =  -4^  and  x  =  0  (see  Fig.  1-3).  An  adjacent,  one¬ 
dimensional  slab  of  liner  material  lies  between  x  =  0  and  x  =  The 
matrix  material  is  characterized  by  a  diffusion  coefficient,  D^,  and  an 
evolution  rate  per  length,  G  e-^,  for  HF  gas.  The  concentration  of  HF  gas 
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in  the  matrix  is  Cj(x,t),  where  t  is  time.  Similarly,  the  HF  concentration 
in  the  liner  is  C2(x,t).  The  liner  is  characterized  by  a  diffusion 
coefficient,  D2,  and  a  rate  of  gas/liner  reaction,  -Hc2(x,t),  which  reduces 
the  value  of  D2,  as  discussed  in  Section  II.  The  initial  values  for  c^  and 
c2  are  zero  for  <  x  <  The  boundary  value  at  x  =  -t1  is  that  each 
incident  molecule  is  lost  to  the  inert  region,  never  to  return.  (The 
opposite  case,  where  no  molecules  are  lost  to  the  inert  region,  is  also 
straightforward  but  is  not  considered  here.)  The  boundary  value  at  x  =  *2 
is  that  every  molecule  that  arrives  at  the  metal  substrate  reacts  with  the 
metal  and  is  immobilized.  The  continuity  condition  at  x  =  0  is  that  the 
fluxes,  Jj  =  -D^3xc^,  and  the  concentrations,  c^,  are  continuous. 

The  barrier  effectiveness,  c,  is  defined  by 

tt> 

f  J2U2,t)  dt 

C  =  1 - 5 - - -  «  100%  ( 18a) 

i  J  J  Ge~gt  dt  dx 
-i,  o 


The  denominator  in  c  is  the  total  number  of  HF  particles  arriving  at  the 
metal  substrate  in  an  unlined  system,  i.e.,  half  the  total  number  of 
evolved  particles,  since  half  of  them  "get  lost"  in  the  inert  material 
(sink).  The  barrier  effectiveness,  t,  varies  between  0%  (when  half  of  the 
evolved  HF  impinges  upon  the  metal  substrate  Just  as  if  no  liner  were 
present)  and  100%  (when  none  of  the  HF  ever  arrives  at  the  metal  sub¬ 
strate).  We  have  to  solve  the  partial  differential  equations  suggested  by 
Fig.  1-3  in  order  to  calculate  c. 


5^4 


The  mathematical  problem  consists  of  partial  differential  equations 
(PDE),  boundary  conditions  (BC),  and  initial  ,'&nditions  (IC): 


PDE 


V%  *  0e"6t 

D?V2c2  -  Hc2  = 


=  atc1 


atC2 


-i1  <  x  <  0 
0  <  x  <  12 


(19a) 

(19b) 


BC 


'D13xc1]o  =  D23x°2]o 


c^Ojt)  =  c2(o,t) 


c^-l^t)  =  0 


c2U2,t)  =  0 


t  >  0 


(20a) 

(20b) 

(20c) 

(20d ) 


c  i  ( x ,  o )  =  0 

-l1  <  x  <  0 

(21a) 

c2(x,o)  =  0 

0  <  x  <  12 

(21b) 

Taking  the  Laplace  transforms  of  Eqs.  (19a)-(19b)  and  using  Eqs.  (21a)  and 
(21b),  we  obtain 


„2-  2- 
V  c,  -  «,c, 

“G 

D^s  +  g) 

(22a) 

V  c2  -  W2C2 

=  0 

(22b) 

“1 

=  /s7d1 

(22c) 

w2 

=  /(s  +  H)/D2 

(22d) 
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Equations  (22a)- (22b)  have  the  solutions 

—  n 

cAs )  =  A,  cosh  w,  x  +  B,  sinh  u.x  -7 - r  (23a) 

11111  s(s  g) 

c2(s)  =  A2  cosh  «i:2x  B2  sinh  u>2x  (23b) 

Applying  Eqs.  (20a)-(20d)  to  Eqs.  (23a)-(23b),  we  obtain 

=  D2w2  coshw2i2sinh(u^l^  +  coshu^t1  sinhu^lp  (24a) 

cd  *  Ai  =  J(f~+  ~^j  (d2u2  coshw2ll2  siHb*0! +  !>.,<*>.,  sinhw^)  (24b) 
-G 

Cd  *  B1  =  ""■■■  j  •  b2< i)2  coshu2i2  (coshu>1l1  -  1)  (24c) 

Q 

cd  *  a2  =  s~(7 ~  g|  ‘  sinhu2t2  (coshw^  -  1)  (24d) 

Cd  *  B2  =  F(TTg7  ‘  Dl“l  coshu)2*2  <coshVr  1}  {24e) 

The  functions,  Eqs.  (23a)-(23b),  with  the  coefficients,  Eqs.  (24a)- 
(24e),  are  the  system  response  functions  on  the  complex  s-plane.  In  order 
to  facilitate  the  following  calculation,  Eqs.  (23a)-(23b)  may  be  written 
more  succinctly  as 


c  As)  =  — - — r 

1  s(s  +  g) 


A.(s)  coshw.x  +  B.(s)  sinhw.x 
_i _ i _ i _ 1  .  1 

Cd(s)  *  6iJ 


(25a) 


where  i  =  1,2,  6^  is  the  Kroneeker  delta  (1  for  i  =  j  and  0  for  i  =  J),  uk 
is  given  by  Eqs.  (22c)-(22d),  Cd(s)  is  given  by  Eq.  (24a),  and 

A ^  =  D2w2  coshw2t2  jinhw^^  +  sinhw2l2  (25b) 

B1  =  ^2^2  coshw2!'2  (coshwilli  -  1) 


(25c) 


A^  =  sinhu^l^  (coshu^A^  -  1) 


(25d) 


B2  =  D1u»1  coshi^tg  (coshw^fc1  -  1)  (25e) 

We  observe  that  Jig  c^(s)  is  finite  for  either  i  =  1  or  2.  (Note  that 

there  are  cancelling  divergent  terns  from  and  the  Kronecker  delta.) 

Hence  there  is  no  pole  in  c, (s)  at  the  origin.  Similarly,  lim  c.(s) 

4  ^  S^^g  1 

diverges  as  (s  ♦  g)  .  Therefore,  there  is  a  simple  pole  at  s  =  -g.  There 
are  also,  potentially,  poles  for  values  of  s  that  cause  Cd(s)  =  0.  We  note 
first  that  s  =  0  and  s  =  -H  cause  Cd(s)  =  0  but  that  Jig  c^(s)  and 
Jin^c^s)  are  nevertheless  finite.  Therefore,  the  roots  s  =  0  and  s  =  -H 
of  Cd(s)  =  0  must  be  rejected,  i.e.,  there  are  no  poles  at  the  origin  (seen 
previously)  or  at  s  =  -H.  All  the  remaining  roots  of  Cd(s)  =  0  are,  in 
general,  simple  poles  on  the  negative,  real  s  axis  (see  Ref.  3)  and  are 
denoted  by  s  =  sn,  n  =  1,2,3,  _ 

There  is,  however,  an  exceptional  condition  on  the  poles  of  c^s)  when 
s  :  sn  :  -g  by  coincidence.  In  this  case,  Jim  c^s)  diverges  as  (s  - 
sn)'2,  i.e.,  there  is  a  second-order  pole  at  s”.  In  this  analysis,  we 
restrict  our  choice  of  g  so  that  this  "coalescence"  of  poles  never  occurs. 
This  simplifies  the  calculation  without  any  significant  loss  of 
generality.  If,  by  any  chance,  the  required  gas  evolution  rate  were  the 
same  as  one  pole  of  the  system  response  function,  second-order  residues 
could  be  calculated.  The  poles  are  summarized  in  Fig.  I I I —1 . 

In  order  to  determine  the  response  functions  c^(s),  it  is  necessary  to 
numerically  determine  the  roots  (sketched  in  Fig.  I I 1-2)  of  Cd(s)  =  0 
(Eq.  (24a),  where 


w2  =  1  for  Is!  >  H 


(26a) 


(26b) 
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Fig.  1 1 1 - 1 .  Pole  Structure  of  the  Response  Function  c.(s).  All  the  poles 
are  simple  if,  as  assumed,  sn  i  -g  for  all  n.  There  are  no 
poles  at  s  =  0  or  s  =  -H. 


III-2.  Roots  of  C^Cs)  =  0  (excluding  s  =  0  and  s  =  -H)  Are  Simple 
Poles  of  tne  Response  Function  c^s).  These  roots  must  be 
found  numerically. 


/H  -  1 3 1 
“2  '  /  D2 


for  | s |  <  H 


(26c) 


When  the  roots  are  known,  c^s)  can  be  expanded  according  to  the  Mittag- 
Leffler  theorem: 


c^s) 


Reslc^s),  -g]  -  Res[c.(s),  si 

- i - +  y  - i - — 

3  *  «  nil  3  "  8n 


(27) 


where  Res[ci(s),  sQ]  refers  to  the  residue  of  ci(s)  at  s  =  sQ. 

Laplace-inverting  Eq.  (27),  we  obtain  the  response  function  in  the 
time  domain  as 

O) 

c ,  { x ,  t)  =  Res[c.(s) ,  -g]  e"gt  *  £  Res[c.(s),s  ]  e^1  (28a) 

1  1  n=1  1  n 


where 


J  r  A.  coshw.x  +  B.  sinhw.x 

Res[c.(s),  -g]  =  =|  (-i - VtTT - -  -  (28b) 


Restc^s)^]  = 


Cd(s) 


X^s) 


11Js  =  -g 


Cd{8) 


s  =  s 


(28c) 


-G 


Xi(s)  =  s^s  ^  ^  [Aa  coshw^x  ♦  B.^  sinhw^x  ♦  6^  •  Cd(s)]  (28d) 


coshw_l_  sinhw.t,  coshw.t,  sinhw^l 


c»  =  \  i - *4 


*ri 


11  .  2  2' 


l1  +  l2 


♦  ( - g - )  sinhu2t2  sinhw^t1 


1  f^lD2“2 

4  5  l  n 


Vi 


Wi 

D2u2 


)  coshu^^  coshu^^ 


(28e) 


When  evaluating  Cqs.  (28a)-(28e),  the  replacements,  Eqs.  (26a)-(26c), 
should  be  performed  separately  for  the  different  |sn|  >  H  and  |sn|  <  H 
regimes.  Using  cosh  ix  =  cos  x  and  sinh  ix  =  1  sin  x,  we  may  keep  track  of 
the  phases  explicitly.  Of  course,  c^(x,t)  must  be  purely  real  and  non¬ 
negative.  In  order  to  calculate  c,  the  spatial  derivative  of  c2(x,t)  must 
be  calculated  [see  Eq.  (18b)].  This  operation  is  performed  on  Eq.  (28a), 
where  the  residues  depend  on  x. 

C.  COMPUTER  CALCULATION  FOR  THE  BARRIER  EFFECTIVENESS,  t. 

OF  FILLED  POLYMERIC  LINERS 

The  equations  from  the  preceding  section  were  programed  on  the 
"DIRAC"  VAX  of  the  Aerospace  Space  Sciences  Laboratory  (SSL) .  Several 
matrix/liner  systems  were  investigated,  as  shown  in  Table  1 1 1- 1 .  Eight 
cases,  defined  in  Table  III-1,  were  investigated,  each  case  consisting  of 
12  possible  combinations  of  parameters. 

The  time-dependent  concentration  profiles  [c^x.t)  and  c2(x,t)]  and 
cumulative  fluxes  (total  evolved,  incident  at  the  inert  sink,  and  incident 
at  the  metal  substrate)  are  shown  explicitly  for  three  of  the  parameter 
combinations  in  Case  I  (Figs.  II 1-3  to  III -8 )  as  well  as  for  three  of  the 
parameter  combinations  in  Case  II  (Figs.  Ill -9  to  III— 14) .  The  barrier 
effectiveness  results  for  the  eight  cases  are  summarized  in  Tables  III —2  to 
III-5.  A  detailed  discussion  of  these  results  follows. 

Figure  II 1-3  shows  the  special  case  in  which  -ly  <  x  <  t2  may  be 
considered  as  one  layer,  because  =  D2  and  H  =  0.  In  this  case,  the 

Sturm-Liouville  theory  is  applicable,  and  the  poles  must  be  in  the  ratio  of 
the  squares  of  the  positive  integers.  Appendix  A  shows  that  the  present 
numerical  calculation  does  give  this  result.  Figure  1 1 1 -3  exhibits  only  a 
very  slight  (i.e.,  imperceptible)  left-right  asymmetry,  because  the  HF  is 
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Table  I 11-1.  Parameters  Used  in  the  Calculation  of  the  Barrier 
Effectiveness,  c,  of  Several  One-Dimensional 
Matrix/Liner  Systems 

1 .  Thicknesses 

matrix:  -i^  s  x  <  0  liner:  0  <  x  < 

2.  HF  generation  rate  in  matrix  =  C  e~ ^ 

[G]  =  particles/(cm  •  Msec) 

[g]  =  Msec-1 

3.  HF  reaction  rate  with  liner  =  -Hc2(x,t) 

[H]  =  Msec"1 
[c2]  =  particles/cm 

4.  In  all  cases,  G  =  10  particles/(cm  •  Msec) 

g_1  =  100  Msec  =  3.17  year 

G*, 

total  HF  s  — —  =  500.0  or  277.0  particles 

5 

depending  on  the  matrix  thickness 

5.  Diffusion  coefficients  matrix:  D1  liner:  D2 

6.  In  all  cases,  H  =  0.0,  10,  100,  and  1000  Msec"1,  in  sequence. 

7.  Cases  investigated 


Case  I : 

D^  =  1.0  cm^/Msec 

S 

=  0.277  cm 

d2  =  1.0,  0.1,  0.01 
(see  Figs.  II 1-3  to  III-8 
and  Table  I I I -2 ) 

l2 

=  0.025 

Case  II: 

D1  =  10.0 

£1 

=  0.277 

D2  =  10,  1.0,  0.1  =  0.025 

(see  Figs.  II 1-9  to  1 11-17 
and  Table  I II 1-2) 
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Table  I 11-1 .  Parameters  Used  in  the  Calculation  of  the  Barrier 
Effectiveness,  c,  of  Several  One- Dimensional 
Matrix/Liner  Systems  (Continued) 


Case 

III: 

D,  =  1.0 

l1  =  0.277 

D2  =  1.0,  0.1,  0.01 

l2  =  0.005 

(see  Table  II 1-3) 

Case 

IV: 

D1  =  10 

l,  =  0.277 

D2  =  10,  1.0,  0.1 

l2  =  0.005 

(see  Table  II 1-3) 

Case 

V: 

D1  =  1.0 

l1  =  0.500 

D2  =  1.0,  0.1,  0.01 

t2  =  0.025 

(see  Table  II 1-4) 

Case 

VI: 

D1  =  10 

t1  =  0.500 

D2  =  10,  1.0,  0.1 

12  =  0.025 

(see  Table  II 1-4) 

Case 

VII: 

D1  =  1.0 

=  0.500 

Dp  =  1.0,  0.1,  0.01 

l2  =  0.005 

(See  Table  III-5) 

Case 

VIII : 

D.|  =  10 

=  0.500 

Dp  10,  1.0,  0.1 

l2  =  0.005 

(see  Table  II 1-5) 
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Fig.  Ill —5 .  HF  Concentration  vs  Position  at  Several  Times.  Since  D^/D2  = 
10,  the  flux  at  the  metal  (right)  is  less  than  that  at  the 
sink  (left).  Matrix:  x  <  0.  Liner:  x  >  0.  Parameters 
above  figure. 
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Fig.  Ill -6 .  Time-Evolution  of  Total  HF,  HF  at  Sink  (left),  and  HF  at  Metal 
(right),  Corresponding  to  Fig.  Ill -5 -  Parameters  above 
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Fig.  Ill  —7 .  HF  Concentration  vs  Position  at  Several  Times.  Since 

100,  the  flux  at  the  metal  (right)  is  less  than  that  at  the 
sink  (left).  Matrix:  x  <  0.  Liner:  x  >  0.  Parameters 
above  figure. 
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Fig.  I II -8.  Time-Evolution  of  Total  HF,  HF  to  Sink  (left),  and  HF  to  Metal 
(right),  Corresponding  to  Fig.  II 1-7.  Parameters  above 
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Fig.  1 1 1  - 1 1 .  HF  Concentration  vs  Position  at  Several  Times.  Since  D.,/D 
10,  the  flux  at  the  metal  (right)  is  less  than  that  at  the 
sink  (left).  Matrix:  x  <  0.  Liner:  x  >  0.  Parameters 
above  figure. 
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Fig.  1 11-12.  Time-Evolution  of  Total  HF,  HF  to  Sink  (left),  and  HF  to 
Metal  (right),  Corresponding  to  Fig.  1 1 1 - 1 1 .  Parameters 
above  figure. 


Table  III-2,  Barrier  Effectiveness,  c  (%) ,  for 
Matrix/Liner  =  0.277/0.025  cm 


D1/D2  H 

(em2/Msec)  (Msec-1) 


0.0 

10.0 

100.0 

1000.0 

1.0/1 .0 

8.3 

8.6 

10.9 

30.0 

1. 0/0.1 

47. *1 

48.5 

57.1 

88.0 

1.0/0.01 

90.0 

91.2 

96.4 

99.97 

10.0/10.0 

8.3 

8.3 

8.5 

10.9 

10.0/1.0 

47.4 

47.5 

48.5 

57.1 

10.0/0.1 

90.0 

90.1 

91.2 

96.4 

Table  111-3. 

Barrier  Effectiveness,  t  (%), 
Matrix/Liner  =  0.277/0.005  cm 

for 

d.|/d2 

2 

(cnr/Msec) 

H 

(Msec-1 ) 

0.0 

10.0 

100.0 

1000.0 

1. 0/1.0 

1.8 

1.8 

1.9 

3.0 

1. 0/0.1 

15.3 

15.4 

16.2 

24.0 

1.0/0.01 

64.3 

64.6 

66.8 

80.7 

10.0/10.0 

1.8 

1.8 

1.8 

1.9 

10.0/1.0 

15.3 

15.3 

15.4 

16.2 

10.0/0.1 

64.3 

64.4 

64.6 

66.8 
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Table  I I 1-4.  Barrier  Effectiveness,  c  (>),  for 
Matrix/Liner  =  0.500/0.025  cm 


Dl'D 2 

H 

p 

(cm  /Msec) 

(Msec-1) 

0.0 

10.0 

100.0 

1000.0 

1.0/1 .0 

4.8 

5.1 

7.6 

27.8 

1.0/0. 1 

33.3 

34.9 

46.9 

86.4 

1.0/0.01 

83.3 

85.4 

94.5 

99.95 

10.0/10.0 

4.8 

4.8 

5.0 

7.6 

10.0/1.0 

33.3 

33.5 

34.9 

46.9 

10.0/0. 1 

83.3 

83.6 

85.4 

94.5 

Table  III-5. 

Barrier  Effectiveness,  t  (It), 
Matrix/Liner  =  0.500/0.005  cm 

for 

H 

p 

(cm  /Msec) 

(Msec-1 ) 

0.0 

10.0 

100.0 

1000.0 

1 .0/1.0 

1 .0 

1.0 

1.1 

2.2 

1. 0/0.1 

9.1 

9.2 

10.2 

18.8 

1.0/0.01 

50.0 

50.4 

53.9 

75.0 

10.0/10.0 

1.0 

1.0 

1.0 

1.1 

10.0/1.0 

9.1 

9.1 

9.2 

10.1 

10.0/0. 1 

50.0 

50.0 

50.4 

53.9 
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generated  only  in  the  region  x  <  0.  The  HF  concentration  rises  to  a 
maximum  before  1  Msec  (278  h  or  11.6  days).  The  maximum  HF  concentration 
occurs  close  to  the  center  of  the  matrix/liner  composite,  where  the  HF 
accumulates  before  it  can  escape  at  the  inert  sink  (at  the  left  of  the 
composite)  or  at  the  metal  substrate  (at  the  right  of  the  composite).  By 
40  Msec  (464  days),  the  maximum  HF  concentration  has  decreased  to  about 
two-thirds  of  its  highest  value.  By  300  Msec  (9.5  years),  the  maximum  HF 
concentration  is  down  to  a  few  percent  of  its  highest  value.  The  diffusion 
relaxation  time  is  0.00924  Msec,  which  is  equal  to  (t^  +  i2)  /Du  .  In 
contrast,  the  characteristic  time  for  HF  evolution  is  g*1  =  100  Msec 
(3.17  years).  The  overall  system  relaxation  time  is  determined  by  the 
longer  of  these  two  times,  i.e.,  the  characteristic  time  for  HF  evolution. 
This  is  clearly  seen  in  Fig.  Ill -3 . 

Figure  II 1-4  shows  the  cumulative  HF  flux  evolved  and  the  fractions  of 
that  cumulative  flux  incident  at  the  inert  sink  and  at  the  metal  substrate. 
The  cumulative  numbers  of  particles  at  the  metal  substrate  and  at  the  inert 
sink  add  up  to  the  total  evolved  number  of  HF  particles  as  time  tends  to 
infinity.  This  is  true  because  H  =  0  and  no  HF  is  bound  up  in  the  system 

for  long  times.  The  cumulative  fluxes  at  the  inert  sink  and  at  the  metal 

substrate  would  be  equal  and  c  would  be  zero  [see  Eq.  (18a)),  except  that 
the  HF  is  generated  only  for  x  <  0.  This  geometric  bias  confers  a  small 
positive  value  on  c,  in  this  case,  8.3V 

Figure  I I I -5  shows  the  case  of  =  1  cm2/Msec,  D2  =  0.1  cm^/Msec,  H  = 
1000  Msec”\  =  0.277  cm,  and  =  0.025  cm.  A  sharp  break  in  the  slopes 
can  be  seen  at  x  =  0,  because  the  fluxes  and  not  the  slopes  are  continuous 
at  that  location.  The  slopes  are  reduced  on  the  metal  substrate  (right) 
side  compared  to  the  inert  sink  (left)  side.  Because  there  is  also  less  D2 

in  the  liner,  the  flux  arriving  at  the  metal  substrate  is  less  than  that 

arriving  at  the  inert  sink  for  two  reasons:  a  smaller  concentration 
gradient  and  a  smaller  diffusion  coefficient.  Appendix  B  indicates  that 
the  poles  are  not  in  the  ratio  of  the  squares  of  the  integers,  because 
D.]  *  D2  and  H  >  0. 
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The  reduction  in  flux  arriving  at  the  wall  is  readily  apparent  in 
Fig.  Ill -6 -  According  to  Eq.  (18a)  or  (18b),  the  barrier  effectiveness,  c, 
is  88%.  Note  that,  in  this  case,  the  numbers  of  HF  particles  having 
arrived  at  the  inert  sink  and  at  the  metal  substrate  at  cumulative  times  do 
not  add  up  to  the  total  number  of  evolved  HF  particles.  This  is  true 
because  H  *  0,  and  some  HF  particles  are  immobilized  in  the  liner. 

The  enhanced  value  of  c  (in  Fig.  I I 1-6  compared  to  Fig.  Ill —4 )  results 
from  a  liner  diffusion  coefficient  (D2)  which  is  less  than  that  of  the 
matrix  (D^),  as  well  as  from  the  reactivity  (H)  of  the  gas  with  the 
liner.  In  Section  II,  it  is  shown  that  H  specifies  a  rate  of  decrease  of 
D2,  which  can  be  interpreted  as  a  progressive  diffusion-channel  clogging. 

A  discussion  follows  of  the  molecular  basis  for  the  existence  of  a 

reduced  D2  in  a  filled  polymer  liner  (see  also  Ref.  12).  At  temperatures 

below  the  polymer  glass  transition  temperature,  T  ,  gas  molecules  must 

© 

diffuse  by  Jumping  between  fixed  free-volume  sites,  which  is  a  relatively 
slow  process.  For  T  >  Tg,  increased  segmental  motion  and  a  higher  free- 
volume  content  allow  faster  diffusion.  Filler  particles  immobilize  polymer 
segments  and  raise  the  effective  Tg  of  the  polymer /filler  composite.  For  a 
given  polymer  at  a  given  operating  temperature,  an  unfilled  system  may  be 
above  its  Tg  (high  D),  while  the  corresponding  filled  system  may  be  below 
its  T^  (low  D) .  Not  all  filled  systems  exhibit  the  drastic  effect  of 
raising  T_  above  the  operating  temperature,  but  any  increase  in  T  results 

O  & 

in  some  reduction  in  D. 

Figure  II 1-7  shows  the  case  in  which  D2  is  down  by  two  orders  of 
magnitude  from  D ^ ,  with  H  =  1000  Msec-1.  In  this  case,  the  slope  disconti¬ 
nuity  in  the  concentration  profile  at  x  =  0  is  sharper,  and  the  concavity 
for  x  >  0  is  more  pronounced.  The  slopes  at  the  metal  substrate  (right 
side)  remain  virtually  zero  at  all  times.  The  success  of  this  liner  in 
keeping  HF  away  from  the  metal  substrate  is  shown  in  Fig.  III-8.  Note 
again  that  the  plateau  values  for  the  cumulative  number  of  particles  at  the 
inert  sink  and  at  the  metal  substrate  as  time  tends  to  infinity  do  not  add 
up  to  the  total  evolved  number  of  HF  particles,  because  H  *  0,  leaving  some 
HF  immobilized  in  the  liner.  Virtually  no  HF  reaches  the  wall  (c  -  99-97%). 
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In  comparison  to  Figs.  II 1-3  to  III —8 1  higher  overall  diffusion 
coefficients  are  shown  in  Figs.  II 1-9  to  III-14.  Figures  II 1-9  and  II 1-10 
are  virtually  identical  to  Figs.  Ill— 3  and  III —4 .  The  reason  is  as 
follows.  The  diffusion  relaxation  time  in  Figs.  III-3  and  III— 4  is 
0.00924  Msec.  The  diffusion  relaxation  time  in  Figs.  II 1-9  and  1 1 1— 10  is 
0.000924  Msec,  being  proportional  to  D”1.  However,  the  overall  system 
relaxation  time  is  dominated  by  the  characteristic  HF  evolution  time,  g" 1  = 
100  Msec.  Hence,  the  change  of  scale  of  D  is  virtually  irrelevant. 

In  contrast,  the  case  in  Figs.  Ill— 1 1  and  III-12  shows  a  significant 
difference  compared  to  Figs.  I II— 5  and  III-6,  even  though  the  ratio  of 
matrix  to  liner  diffusion  coefficients  is  the  same.  The  concavity  of  the 
concentration  profile  within  the  liner  is  reduced,  while  the  slopes  at  the 
metal  substrate  are  increased.  The  barrier  effectiveness,  c,  is  reduced  to 
57. U  in  Fig.  II 1-12  (it  had  been  88%  in  Fig.  II 1-6) .  Figures  111-13  and 
II I— 14  show  the  concentration  profiles  and  cumulative  fluxes  for  - 

100  (compare  to  Figs.  I I 1-7  and  II 1—8) .  There  is  still  a  fair  amount  of 
concavity  in  the  liner  concentration  profiles.  Only  slightly  more  HF 
arrives  at  the  metal  substrate  (c  =  96.4 %  compared  to  99. 97^  in 
Fig.  1 1 1 -8 ) . 

The  barrier  effectiveness  results  for  all  the  cases  listed  in 
Table  1 1 1 - 1  are  presented  in  Tables  1 1 1-2  to  III-5.  The  liner  thicknesses 
are  either  0.025  cm  (about  10  mils)  or  0.005  cm  (about  2  mils). 

In  the  first  row  (the  horizontal  listing)  of  Tables  III -2  and  III —4 , 
there  is  a  substantial  increase  of  c  with  increasing  H  at  D^/D2  =  1. 
However,  the  maximum  ;  is  only  about  30% .  In  all  the  other  rows  of 
Tables  III -2  to  I I 1—5 ,  there  is  only  a  modest  increase  of  c  with  increasing 
H.  In  contrast,  there  is  always  a  large  c  increase  in  each  column  (the 
vertical  listing)  in  Tables  1 1 1-2  to  III-5,  where  D^/D2  increases  at 
constant  H.  This  shows  that,  of  the  two  effects  caused  by  filler  particles 
(i.e.,  reduced  D2  and  increased  H),  the  D2  reduction  is  far  more  important 
in  increasing  c.  This  result  indicates  that  it  is  useful  to  define  a 
logarithmic  diffusion-coefficient  ratio  as 
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R  =  logCD^Dg) 


(29) 


For  matrix  and  liner  thicknesses  of  0.277  and  0.025  cm,  respectively, 
and  R  =  2,  Tables  111-2  to  III -5  show  that  barrier  effectiveness  in  excess 
of  90>  can  be  achieved.  In  contrast,  matrix/liner  thicknesses  of  0.277/ 
0.005  cm  result  In  ;  >  65J  for  R  =  2.  For  matrix/liner  thicknesses  of 
0.500/0.025  and  0.500/0.005,  the  c  values  are  at  least  83J/50 %  for  R  =2. 

Figure  II 1-15  shows  c  plotted  as  a  continuous  function  of  D2  (1-2  is  a 
parameter).  Most  of  the  increase  in  ;  occurs  for  D2  <  0.1  cm2 /Msec,  corre¬ 
sponding  in  this  case  to  R  >  1.  Finally,  in  Fig.  III-16,  c  is  plotted  as  a 
continuous  function  of  (D2  is  a  parameter).  Larger  c  values  are 
achieved  for  thinner  liners  when  R  =  2  (top  curve)  compared  to  R  =  1  (lower 
curve) . 

D.  COMPUTER  CALCULATIONS  OF  THE  BARRIER  EFFECTIVENESS,  c.  FOR  A 

REALISTIC  CASE  OF  A  FILLED  POLYMER  CONTAINING  DEFECTS 

In  Section  II,  single-slab,  diffusion-profile  measurements  at  two  HF 
exposure  times  led  to  estimates  of  the  diffusion  coefficients  of  an  un¬ 
filled  and  of  a  filled  (AI2O3)  model  polymer  (Uralane  5753).  The  diffusion 
coefficient,  D2,  of  the  unfilled  system  was  estimated  to  be  at  least 
0.8  cm  /Msec  and  of  the  filled  system,  at  least  0.1  cnr/Msec.  In  either 
case,  H  was  at  least  13  Msec-^ .  Until  measurements  at  additional  times  are 

available,  we  provisionally  assign  D2  =  0.8  cm2/Msec  to  the  unfilled 

P  -  1 

system,  =  0.1  cnr/Msec  to  the  filled  system,  and  H  =  13  Msec  in  either 

2 

case.  Further  assuming  that  =  1  cm  /Msec,  which  is  coninon  for  gas 
diffusion  in  typical  polymers  (see  Section  II),  we  can  calculate  the 
barrier  effectiveness,  c,  of  the  model  polymer  in  both  its  unfilled  and 
filled  states.  Of  course  there  are  other  mechanical,  thermal,  and  adhesion 
properties  that  have  to  be  evaluated  in  the  selection  of  an  actual  liner. 
This  calculation  merely  illustrates  what  could  be  achieved  in  an  actual 
system. 
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Fig.  I I I - 1 5 .  Barrier  Effectiveness,  c,  vs  the  Liner  Diffusion  Coefficient, 
D^,  for  Two  Values  of  Liner  Thickness.  Parameters  above 
figure . 


Fig.  1 1 1  —  1 6 .  Barrier  Effectiveness,  c,  vs  the  Liner  Thickness,  t„>  f°r  T”0 
Values  of  the  Liner  Diffusion  Coefficient,  Dg.  Parameters 
above  figure. 


In  Figs.  II 1—17  and  I II— 18,  we  see  that  the  barrier  effectiveness  of 
the  unfilled  system  is  very  small  (c  =  10.6}),  as  expected  for  R  = 
log(Dj/D2)  =  log( 1/0.8)  =  0.097.  Because  of  the  weak  H  dependence  noted 
above,  the  system  may  be  considered  to  exist  at  D2  =  0.8  on  the  top  curve 
in  Fig.  II I- 15 .  In  Figs.  II I- 19  and  III-20,  we  see  that  c  for  the  filled 
system  is  much  larger  (c  =  48.9}),  because  R  =  1 .  The  system  may  be 
considered  to  exist  at  D2  =  0.1  on  the  top  curve  in  Fig.  Ill— 15  and  at  = 
0.025  (extreme  right)  on  the  lower  curve  in  Fig.  I II— 16.  These  initial 
results  possibly  indicate  that  actual  liner  systems  can  be  designed  with 
R  ♦  2  and  c  ♦  90}. 

Objections  have  been  made  that  the  experimental,  single-slab  HF 
concentration  profiles  of  the  Al20g-filled  polymer  in  Section  II  do  not 
look  exactly  like  the  computed  curves  in  Fig.  II 1-21 .  Specifically,  the 
experimental  data  show  a  low-level  plateau  (for  x  >  0.3  cm),  which  is  not 
evident  in  the  computed  model  response.  In  contrast,  the  experimental 
response  of  the  unfilled  system  does  match  the  computed  curves.  The  low 
level  of  additional  diffusion  is  attributed  to  some  unspecified  defect  in 
the  filled  system.  One  possible,  but  unverified,  explanation  is  that 
cracks  have  developed  between  filler  and  polymer,  allowing  additional 
diffusion.  It  has  been  speculated  that  this  could  be  a  "fatal"  drawback  to 
the  use  of  filled  polymers  as  liners.  He  show  in  the  following  discussion, 
however,  that  the  observed  effect  is  consistent  with  only  a  small  reduction 
in  the  barrier  effectiveness,  c* 

The  effect  of  the  defects  can  be  modeled  by  assuming  that  some  frac¬ 
tion  of  the  diffusing  species  finds  faster  diffusion  pathways,  which  may  be 
traversed  in  parallel  with  diffusion  through  the  polymer.  The  single-slab 
response  of  such  a  faster  pathway  is  shown  in  Fig.  I I 1-22,  where  D  =  1  and 
H  =  0  in  contrast  to  Fig.  II 1-2 1  (D  =  0.1  ard  H  =  13).  It  is  essential 
that  H  =  0  in  such  a  pathway,  so  that  the  slope  is  zero  at  long  times,  as 
observed  experimentally. 

Considering  parallel  pathways  in  which  90}  of  the  particles  diffuse 
through  the  polymer  and  10}  through  the  defect  channels,  we  compute  the 
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Fig.  1 1 1 - 1 7 -  HF  Concentration  vs  Position  at  Several  Times,  Corresponding 
to  an  Unfilled  Model -Polymer  Liner  as  Measured  in  Section 
II.  Parameters  above  figure. 


Fig.  1 1 1 - 1 8 .  Time-Evolution  of  Total  HF,  HF  to  Sink  (left),  and  HF  to 
Metal  (right),  Corresponding  to  Fig.  III-17.  c  =  10. 6%. 
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III —2 1 .  Single-Slab  Profiles  of  HF  Concentration  Which 
Represent  the  A^O^-Filled  Uralane  Data  in 
Fig.  1 1- 13b.  For  x  >  0.3  cm,  the  experimental 
data  show  a  low-level  plateau  not  evident  in  this 
calculation. 


D  H  t 


Fig.  1 11-22.  Single-Slab  Results  for  the  Parameters  Given  Above  Figure. 

Since  H  =  0,  the  concentration  approaches  unity  everywhere 
for  long  time. 
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single-slab  response  (see  the  equations  in  Section  II)  in  Pig.  III-23. 

This  result  adequately  represents  the  experimental  data  in  Fig.  II-13b. 
Actually,  the  representation  is  believed  to  be  conservative,  because  the 
low-level  plateau  extends  farther  in  the  computed  response. 

This  same  90J/10J  ratio  of  diffusing  particles  in  the  (D2  =  0.1,  H  = 
13)/(D2  =  1,  H  =  0)  channels  can  be  used  in  the  double-layer  calculations 
Just  described.  The  results  are  shown  in  Figs.  II 1-24  and  III -25 .  There 
is  only  a  small  difference  from  the  results  shown  in  Figs.  I II- 19  and 
III-20.  In  fact,  c  =  HH.8%  in  the  defective  system,  compared  to  48. 9 %  in 
the  perfect,  filled  system.  This  result  is  also  intuitively  clear,  because 
the  low-level  plateau  has  essentially  zero  slope  and  the  flux  is  propor¬ 
tional  to  the  slope.  Of  course,  it  is  possible  to  imagine  the  formation  of 
giant  cracks  in  a  filled  polymer  exposed  to  HF  and  the  drastic  reduction  of 
its  barrier  effectiveness.  However,  this  view  is  contradicted  by  the 
available  experimental  evidence.  The  observed  defects  (Section  II)  are  not 
a  catastrophic  drawback  to  the  use  of  filled  polymers  as  barrier  liners. 

E.  SUMMARY 

The  diffusion-barrier  properties  of  several  matrix/liner  systems  were 
investigated  as  a  function  of  thickness  (matrix  =  0.277  or  0.500  cm,  and 
liner  =  0.025  or  0.005  cm);  diffusion  coefficient  (matrix  =  10  or 

y  P 

1.0  cnr/Msec,  and  liner  D2  =  10,  1,  0.1,  or  0.01  cnr/Msec) ;  and  reactivity 
(H  =  0,  10,  100,  or  1000  Msec"1).  The  properties  were  investigated  in  a 
system  in  which  HF  gas  is  generated  with  a  characteristic  time  of  100  Msec 
(3.17  years).  All  parametric  combinations  are  defined  in  Table  1 1 1 - 1 . 

Some  examples  of  the  time-dependent  concentration  profiles  and 
cumulative  fluxes  are  shown  in  Figs.  Ill -3  to  II I- 14 .  The  results  for  the 
barrier  effectiveness,  c,  of  all  the  parametric  combinations  are  given  in 
Tables  III -2  to  III-5,  where  c  is  defined  in  Eqs.  ( 1 8a ) - ( 1 8b )  (see  also 
Fig.  1-3).  The  barrier  effectiveness,  c,  is  zero  when  half  of  the  evolved 
HF  particles  impinge  on  the  metal  substrate  and  the  other  half  escapes  to 
the  inert  sink.  The  barrier  effectiveness,  c,  is  100J  when  none  of  the 
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X  (cm) 


Fig.  Ill ~23 -  Parallel-Pathways  Model  Result  (single  slab)  for  a  Filled 

Polymer  (D  =  0.1,  H  =  13)  with  Defects  (D  =  1,  H  =  0).  The 
polymer/defect  conductance  ratio  is  90%/ 10% .  This  result 
adequately  represents  the  experimental  data  in  Fig.  I I - 1 3b . 
The  defects  may  be  provisionally  thought  of  as  cracks. 


X  (cm) 


Fig.  111-24.  Parallel -Pathways  Result  (double  slab)  for  a  Filled  Polymer 
Liner  (D  =  0.1,  H  =  13)  with  Defects  (D  =  1,  H  =  0)  in  the 
Region  x  >  0.  The  polymer /defect  conductance  ratio  is 
90/&/105L  The  matrix  D,  =  1  in  the  region  x  <  0.  This  result 
should  be  compared  to  Fig.  1 11-19. 


EVOLUTION 


Fig.  Ill -25 .  Time-Evolution  of  Total  HF,  HF  at  Sink  (left),  and  HF  at 
Metal  (right),  Corresponding  to  Fig.  Ill -24  (defective 
liner).  Compared  to  Fig.  111-20  (no  defects),  we  see  that 
the  barrier  effectiveness  is  only  reduced  from  48 . 9%  to 
44 . 8% .  These  defects  may  be  provisionally  thought  of  as 
cracks. 
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evolved  HF  impinges  on  the  metal  substrate,  instead  either  escaping  to  the 
inert  sink  or  being  inmobilized  in  the  liner. 

Only  a  relatively  weak  dependence  of  5  on  the  reactivity,  H,  was 
noted.  However,  filler  particles  in  real  polymer  systems  also  serve  to 
reduce  the  liner  diffusion  coefficient,  and  this  reduction  can  have  a 
dramatic  effect  on  t.  For  the  case  of  the  matrix  =  0.277  cm  and  the 
liner  i ^  =  0.025  cm,  c  assumes  values  in  excess  of  90%  for  R  =  log(D1/D2)  = 
2  and  values  in  excess  of  M7%  for  R  =  1.  For  a  given  parametric 
combination,  t  increases  monotonieally  with  R  (the  logarithmic  diffusion- 
coefficient  ratio) .  The  most  beneficial  computed  effects  were  obtained  for 
R  =  2,  which  should  be  the  goal  of  future  experimental  liner  development. 

The  detection  of  defects  (slightly  nonclassical ,  experimental  concen¬ 
tration  profiles)  in  filled  polymers  exposed  to  HF  raised  the  possibility 
that  the  beneficial  effects  of  the  filler  could  be  "short-circuited"  by  the 
introduction  of  faster,  parallel  diffusion  pathways.  The  precise  nature  of 
these  defects  has  not  been  determined,  but  they  may  be  provisionally 
thought  of  as  cracks  between  the  filler  particles  and  the  polymer.  It  was 
shown,  however,  that  (1)  the  measured  effect  of  the  defects  may  be  success¬ 
fully  modeled  by  assuming  that  10%  of  the  HF  travels  in  the  defect  pathways 
and  (2)  the  barrier  effectiveness  of  such  a  defective  liner  (based  on 
measurement  in  Section  II)  is  only  reduced  from  48. 9%  to  44.8%  by  these 
defects.  At  this  time,  there  is  no  observable  fatal  disadvantage  to  the 
use  of  a  filled  polymer  as  a  barrier  to  HF. 


IV.  EFFECT  OF  FILLER  TYPE  ON  HF  DIFFUSION 


A.  BACKGROUND 

In  the  previous  section,  we  investigated  the  effectiveness  of  unfilled 
and  alumina-filled  polymer  systems  as  barrier  liners  in  reducing  the  amount 
of  HF  gas  impinging  on  a  metal  substrate.  In  this  section,  freshly 
prepared  Uralane  5753  was  used  alone  and  in  combination  with  magnesium 
silicate,  strontium  chromate,  and  carbon  black.  The  transport  coefficients 
of  these  liner  materials  were  determined  so  that  multislab  calculations 
could  reveal  the  extent  to  which  the  barrier  reduces  HF  flux  at  the  metal 
substrate  (see  Fig.  1-3).  One -dimensional  diffusion  experiments  were  used, 
as  in  Section  II,  to  determine  both  the  basic  diffusion  coefficient,  D,  and 
a  second  parameter,  H,  specifying  the  degree  of  chemical  reactivity  between 
the  gas  and  the  liner  material. 

In  Section  II,  some  experiments  (of  21  h  and  12  day  durations)  on  a 
model  polymer  with  or  without  alumina  particles  (40  wt%  of  10  ±  2  urn 
particles)  were  found  to  be  consistent  with  H  -  13  Msec-1  for  either  filled 

or  unfilled  polymer:  D  '  0.8  cm  2 /Msec  for  unfilled  polymer  and  D  -  0.1 

2  A 

cm  /Msec  for  filled  polymer.  (One  roegasecond  is  10  sec  or  278  h.)  In 

Section  III,  it  was  shown  that  barrier  effectiveness  of  the  unfilled  system 

was  only  10.6%,  while  that  of  the  filled  system  was  44. 8£.  The  barrier 

effectiveness  is  the  percent  reduction  in  the  number  of  HF  particles 

impinging  on  the  wall  due  to  the  presence  of  the  liner. 

The  enhanced  value  of  the  barrier  effectiveness  results  from  a  liner 
diffusion  coefficient  (D2)  which  is  less  than  that  of  the  matrix  ( D ^ ) ,  as 
well  as  from  the  reactivity  (H)  of  the  gas  with  the  liner.  In  Section  II, 
it  is  shown  that  H  specifies  a  rate  of  decrease  of  D2,  which  can  be 
interpreted  as  a  progressive  diffusion-channel  clogging.  We  suggest  that 
the  molecular  basis  for  the  existence  of  a  reduced  D2  in  our  filled 
polymers  is  the  following  (Refs.  12-16):  At  temperatures  below  the  polymer 
glass  transition  temperature,  T  ,  gas  molecules  must  diffuse  by  Jumping 
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r 
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between  fixed  free-volume  sites,  which  is  a  relatively  slow  process.  For  T 

>  Tg,  increased  segmental  motion  and  a  higher  free-volume  content  allow 

faster  diffusion.  Filler  particles  immobilize  polymer  segments  and  raise 

the  effective  Tg  of  the  polymer/filler  composite.  For  a  given  polymer  at  a 

given  operating  temperature,  an  unfilled  system  may  be  above  its  T_  (high 

© 

D),  while  the  corresponding  filled  system  may  be  below  its  T'  (low  D).  Not 

S 

all  filled  systems  may  exhibit  the  drastic  effect  of  raising  T  above  the 

O 

operating  temperature,  but  any  increase  in  Tg  will  result  in  some  reduction 

in  D.  Work  is  in  progress  to  determine  the  T  for  the  filled  and  unfilled 

D 

Uralane  under  consideration  here. 

In  this  section,  the  relative  effects  of  filling  the  polymer  with 
magnesium  silicate,  strontium  chromate,  and  carbon  black  will  be  assessed 
by  calculating  the  barrier  effectiveness,  c,  which  is  defined  in 

Section  III  as  the  percent  reduction  in  the  number  of  HF  particles  imping¬ 

ing  on  the  metal  substrate  due  to  the  presence  of  the  liner. 

B.  DIFFUSION-FRONT  TRAJECTORIES 

Consider  a  one-dimensional  slab  of  material  between  positions  x  =  0 

and  x  =  1.  The  slab  is  exposed  to  a  constant  gas  concentration,  c  ,  at  x  = 

2 

0.  Diffusion  of  gas  into  the  slab  proceeds  according  to  D3xc  -  He  =  3tc, 
where  D  is  the  diffusion  coefficient  and  H  is  the  rate  of  reaction  between 
the  gas  and  the  solid  (see  Section  II).  The  concentration,  c(x,t),  of  gas 
in  the  slab  is  given  by 
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Figure  I V- 1  is  an  example  of  this  solution  for  the  case  in  which  D  = 

P  1 

0.5  cm  /Msec,  H  =  13  Msec  ,  and  l  =  1  cm  at  the  successive  times  0.001, 
0.004,  0.007,  0.010,  0.025,  0.040,  0.070,  0.100,  and  0.500  Msec. 

For  the  case  in  which  H  =  0  and  t  ♦  «,  this  solution  may  also  be 
written  as  (see  Ref.  19): 


where 


SiiUi  ,  ,  .  erf(_i_) 

o  2/Dt 


(3D 


erf  (z)  =  ^7- 


z  2 
f  e_T1  dn 


(32) 


Consider  c(x',t'),  where  c(x',t')/c0  =  0.1573.  In  that  case,  the  argument 
of  erf  (z)  is  unity,  and 


x’  =  2  /Dt’  (33) 

at  c(x,t)/cQ  =  0.1573 
for  H  =  0  and 

In  Section  II,  it  was  hypothesized  that  this  relation  could  be  generalized 
to 


x’  =  Jp-  erf  (Hf )  (34a) 

Ht'  <<  1  (34b) 

t'  ♦  «  (34c) 

for  l  -►  ®.  In  Section  II,  x'  was  determined  from  c(x',t')/c0  =  0.17  for 
convenience. 

The  concentration  profiles  in  Fig.  IV-1  are  all  very  small,  at  x  =  l, 
so,  effectively,  i  -  ®.  These  profiles  are  replotted  in  Fig.  I V -2 ,  and  the 


x'  •  2  /Dt' 


87 


1.02 


0.5  Msec 
0.1  Msec 
0.07  Msec 
0.04  Msec 
0.025  Msec 
0.01  Msec 
0.007  Msec 
0.004  Msec 
0.001  Msec 


0.70 


0.80 


0.90 


Fig.  I V- 1 .  Calculated  Concentration  Profiles  for  the  Diffusion  of  Gas 
into  a  Solid,  Including  the  Reactivity  Between  the  Gas  and 
Solid.  Parameters  given  above  figure. 


Fig.  IV-2.  Same  Concentration  Profiles  as  in  Fig.  IV-1  Replotted  on  a 
Smaller  Scale 
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successive  (x',t')  data  points  at  which  c(x',t')/c0  =  0.1573  are  plotted  in 
Fig.  I V-3-  There  is  excellent  agreement  with  Eqs.  (34a)-(3Mc) ,  because  the 
curve  shape  is  like  the  erf(z)  function,  the  slope  (1.M25)  at  the  origin 
gives  D  =  0.51  cm^/Msec  [Eq.  (3^b)],  and  the  plateau  (0.36  cm)  gives  H  = 

12  Msec-1  [Eq.  (3*<c)].  This  "slope/plateau”  approach  is  the  method  for 
determining  the  fundamental  gas-transport  quantities  D  and  H  from  the 
experimental  concentration  profiles. 

It  is  beyond  the  scope  of  this  investigation  to  measure  an  entire 
trajectory  like  that  of  Fig.  I V-3 -  To  provide  an  approximate  evaluation, 
the  following  approach  was  adopted  in  Section  II.  In  the  cases  of  unfilled 
and  A^O^-filled  Uralane,  it  was  noted  that  there  was  no  further  advance  of 
the  diffusion  front  at  288  h  compared  to  21  h  of  HF  exposure.  If,  on  the 
other  hand,  H  had  been  2ero,  the  diffusion  front  would  have  advanced  by  a 
factor  of  3-7,  according  to  Eq.  (33).  If  the  generalization,  Eq.  (3^a),  is 
valid,  then  this  immobilization  of  the  diffusion  front  should  occur  after  a 
characteristic  time,  H-1.  If  this  time  is  21  h  or  less,  then  H  is 

13  Msec*1  or  more.  It  is  considered  conservative  to  set  H  at  the  lower  end 
of  this  range  without  jeopardizing  the  subsequent  calculations  of  barrier 
effectiveness,  which  are  much  less  sensitive  to  H  than  to  D  (see 

Section  II).  Therefore,  in  all  cases  analyzed,  we  set  H  to  13  Msec-1  as  a 
reasonable  estimate  until  an  entire  trajectory  can  be  measured. 

Using  this  value  of  H,  the  value  of  D  for  a  gas/polymer  system  may  be 
obtained  from  the  long-time  (stationary)  concentration  profiles  as  time 
tends  to  infinity,  as  follows.  From  Eq.  (30),  we  obtain  the 
slope,  oQ,  where 

|co|  =  lim  8xc|x=0  tanh  t)  (35) 

Measurement  of  oQ  and  knowledge  of  H  and  l  allow  D  to  be  calculated. 
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C.  EXPERIMENTAL  METHOD 


Polyethylene  cups  of  unfilled  and  filled  Uralane  5753  were  freshly 
prepared  for  one-dimensional  HF  diffusion  experiments,  as  described  in 
Section  II.  The  fillers  were  magnesium  silicate  (2Mg0’3SiC>2’xH20,  Lot  No. 

CE  311,  from  Spectrum  Chemical,  with  initially  20  wt$  MgO  and  45  wt£  Si02); 
strontium  chromate  [SrCrOjj,  Lot  No.  D13C,  from  Alfa  Products /Morton  Thiokol 
(98% ) ] ;  and  carbon  black  (National  Spectroscopic  Powder,  Lot  No.  G66,  Grade 
SP2) .  The  fillers  were  baked  at  120°C  for  24  h  before  their  incorporation 
(40  wt>)  into  the  room-temperature  cured  Uralane.  The  Uralane  and  filler 
mixture  was  stirred  by  hand  until  uniform  to  the  eye.  Since  the  typical 
XPS  beam  dimension  is  about  0.5  mm,  while  the  average  interparticle  spacing 

ft  o 

is  30  ym  (or  less),  the  XPS  averages  over  an  area  containing  17  =  289  or 

more  particles.  This  being  the  case,  variations  in  interparticle  spacing 
due  to  mixing  are  not  expected  to  be  observable  in  the  XPS  experiment.  The 
generally  monotonic  results  for  fluorine  concentration  vs  depth  confirm 
this  expectation.  Curing  continued  for  2  weeks  before  the  HF  exposure 
began.  The  Uralane,  part  A,  was  from  batch  GLL-017,  while  part  B  was  from 
batch  HM-056.  Both  batches  were  dated  January  1988. 

The  Le  Mont  particle  analyzer  was  used  to  determine  the  filler 
particle-size  distributions  for  the  Mg0‘Si02  and  the  SrCrO^.  There  was 
insufficient  contrast  to  the  carbonaceous  substrate  to  allow  this  scanning 
electron  microscopy  (SEM)-based  technique  to  be  used  for  the  carbon 
black.  A  different  technique  might  allow  results  for  the  same  grade  of 
carbon  black  to  be  determined  in  the  future.  The  complete  analysis  for 
MgO-SiC^  and  SrCrO^  is  given  in  Appendix  C.  Typical  particle  lengths  are 
4.6  yin  ±  54 %  for  MgO-SiC^  and  5.2  ym  ±  61J  for  SrCrO^. 


If  a  40  wt.%  sample  is  considered  (conservatively)  to  have  4  volj  of  5  pm 
particles,  then  it  has  42/3  =  2.5  area  percent  of  particles.  Each  25  pm2 
particle  therefore  exists  on  an  average  area  of  1000  pm2.  The  average 
interparticle  distance  is  therefore  ^1000,  or  approximately  30  pm 
(conservatively).  If  the  volume  percent  is  larger,  then  the  interparticle 
spacing  is  smaller. 


The  samples  were  exposed  to  HF  for  9  days  in  a  large  polyethylene 
chamber.  The  presence  of  the  polyethylene  cups  around  the  samples  forced 
all  the  HF  diffusion  to  proceed  from  one  face  of  each  sample.  After 
exposure,  the  samples  were  sectioned  for  x-ray  photoelectron  spectroscopy 
(XPS)  analysis.  The  fluorine  concentration  normalized  by  the  carbon 
concentration  was  recorded. 

D.  EXPERIMENTAL  RESULTS 

Figure  IV-4  shows  the  carbon-normalized  fluorine  concentration  vs 

penetration  depth  for  the  unfilled  Uralane  after  its  9-day  exposure  to 

HF.  The  solid  line  is  a  least-squares  curve  fit  [c(x)  =  a  -  oQ  x]  to  the 

data,  excluding  the  points  at  x  =  1.9  and  3.8,  which  originate  from  regions 

considered  to  be  smeared  during  sample  preparation.  The  theoretical 

response  function  over  the  curve-fit  region  is  not  linear;  the  curve  fit  is 

merely  a  way  to  objectively  estimate  the  slope  at  x  =  0.  It  is  found  that 

(a,  o  )  =  (1.720,  0,2456)  for  c(x)  in  percent  and  that  (a,  o  )  =  (1.0, 
o  o 

0.11128  ran-1)  for  normalized  c(x)/cQ.  Since  the  sample  is  about  25  urn  deep, 
we  may  assume  t  •*  =>.  Hence,  o  =  1.428  cm'1  =  /H?D,  from  Eq.  (35). 

Assuming  H  =  13  Msec-1,  as  just  noted,  D  =  6.4  cnr/Msec.  This  result  for 
freshly  prepared  Uralane  is  considerably  larger  than  that  obtained  for 
several-year-old  Uralane  (D  =  0.8  cm  /Msec)  in  Section  II.  Work  is  in 
progress  for  the  determination  of  Tg  for  aged  samples.  As  noted  in  Section 
IV. A,  there  is  expected  to  be  a  correlation  between  T„  and  the  diffusion 

D 

coefficient. 

Figure  IV-5  shows  the  concentration  profile  for  the  magnesium- 
silicate-filled  Uralane  after  its  9-day  HF  exposure.  Since  the  concentra¬ 
tion  rapidly  decreases  as  x  increases,  the  concentration  drops  below  the 
XPS  detection  limit.  Thus,  the  last  two  data  points  in  Fig.  IV-5  show  zero 
concentration. 

The  solid  line  in  Fig.  IV-5  is  a  least-squares  curve  fit  to  the  data, 
excluding  the  first,  seventh,  and  eighth  data  points.  The  first  data  point 
Is  obviously  smeared,  while  the  last  two  data  points  are  not  on  the  approx- 
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Fig.  IV-5.  Experimental  Fluorine  Concentration  Profile  for  MgO-SiC^- 
Filled  Uralane  After  9  Days  of  HF  Exposure  " 
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Fig.  IV-6.  Experimental  Fluorine  Concentration  Profile  for  SrCrO^-Filled 
Uralane  After  9  Days  of  HF  Exposure 
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imately  linear  portion  of  the  profile.  It  is  found  that  (a,  aQ)  =  (23.05, 
6.910)  for  c  in  percent  and  that  (a,  oQ)  =  (1.0,  0.2998  mm-1)  for  normal¬ 
ized  c.  Using  oQ  =  /H/D  and  H  =  13  Msec*1,  we  find  D  =  1.4  cm2/Msec. 

Thus,  these  4.6  ym  (t  54?)  magnesium  silicate  particles  reduce  D  by  a 
factor  of  4.5. 

Figure  IV-6  shows  the  concentration  profile  for  the  strontium- 
chromate-filled  Uralane  after  its  9  days  of  HF  exposure:  The  solid  line  is 
a  least-squares  curve  fit  to  the  data,  excluding  the  two  data  points  not  on 
the  linear  portion  of  the  profile.  It  is  found  that  (a,  oQ  )  =  (5.948, 
2.680)  for  c  in  percent  and  that  (a,  oQ)  =  (1,  0.4506  mm*1)  for  normalized 
c.  As  before,  it  is  found  that  D  =  0.64  cm2 /Msec.  Thus,  these  5.2  ym  (± 
60?)  strontium  chromate  particles  reduce  D  by  a  factor  of  10. 

Figure  IV-7  shows  the  concentration  profile  for  the  carbon-black¬ 
filled  Uralane  after  its  9-day  HF  exposure.  Curve  fitting  reveals  that 
(a,  oq)  =  (1.522,  0.7213)  for  c  in  percent  and  that  (a,  oQ)  =  (1.0, 

0.4739  mm*1)  for  normalized  c.  Thus  D  =  0.58  cm2/Msec,  which  is  an  11 -fold 
decrease  from  the  unfilled  condition. 

In  Table  IV- 1 ,  the  preceding  results  are  summarized.  In  addition,  the 
barrier  effectiveness,  {,  of  a  polymer  liner  with  these  D  and  H  values  is 
calculated.  The  barrier  effectiveness,  c,  is  the  percent  reduction  in  the 
number  of  HF  particles  impinging  on  the  metal  substrate  due  to  the  presence 
of  the  liner  and  is  rigorously  defined  in  Section  II.  In  these 
calculations,  the  matrix  thickness  is  0.277  cm,  and  the  liner  thickness  is 
0.025  cm. 

E.  SUMMARY 

We  investigated  the  effectiveness  of  unfilled  and  filled  polymers  as 
barrier  liners  to  reduce  the  amount  of  HF  gas  that  evolves  from  a  matrix 
and  impinges  upon  a  metal  substrate  (see  Fig.  1-3).  The  polyurethane, 
Uralane  5753,  was  used  alone  and  in  combination  with  magnesium  silicate, 
strontium  chromate,  and  carbon  black.  Typical  particle  lengths  were  4.6  ym 
(for  magnesium  silicate)  and  5.2  ym  (for  strontium  chromate).  The  carbon 
black  size  distribution  was  not  determined  because  of  lack  of  contrast 
(relative  to  the  substrate)  in  the  SEM.  After  a  9-day  exposure  to  HF  gas, 
the  fluorine  concentration  vs  depth  was  determined  by  XPS.  Assuming  that 
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the  reactivity,  H,  was  always  13  Msec-1,  as  in  Section  II,  we  determined 
that  the  diffusion  coefficients  for  unfilled,  magnesium-silicate-filled, 
strontium-chromate-filled,  and  carbon-black-filled  Uralane  were  6.4,  1.4, 
0.64,  and  0.58,  respectively. 

The  barrier  effectiveness,  c,  of  a  barrier  liner  has  been  previously 
defined  (Section  III)  and  is  the  percent  reduction  in  the  number  of  HF 
particles  impinging  on  the  motor  wall  due  to  the  presence  of  the  barrier 
liner.  We  assumed  that  the  matrix  has  D  =  6.4  cm^/Msee  and  evolves  HF  at  a 
rate  Ge-^,  where  g  =  0.01  Msec-1  and  G  =  10  particles/(  cm-Msec) .  We 
determined  that  barrier  effectivenesses  are  c  =  8.3,  29.4,  ‘'7.7,  and  50.1%, 
for  the  unfilled,  magnesium-silicate-filled,  strontium-chromate-filled,  and 
carbon-black-filled  Uralane,  respectively.  In  this  calculation,  a  matrix 
thickness  of  0.277  cm  was  used,  and  the  liner  thickness  was  0.025  cm. 
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APPENDIX  A. 


RELATION  OF  COMPLEX  ANALYSIS  TO  STURM-LIOUVILLE  THEORY 

This  appendix  presents  the  poles  of  the  response  function  for  a 
double-layer  system  with  the  given  parameters.  Note  that  in  this  case,  D 
=  D2  and  H  =  0  (see  Figs.  II I— 3  and  111-4),  so  that  the  Sturm-Liouville 
theory  is  applicable  and  the  poles  must  be  in  the  ratio  of  the  squares  of 
the  integers.  These  poles  were  found  numerically  and  are  in  the  proper 
ratio  to  at  least  seven  significant  digits  in  this  double-precision 
calculation. 

Case  I  (Table  I II -1 ) :  D1  =  1.0  cm2 /Msec 

D2  =  1.0 

H  =  0.0  Msec-1 

_J  SN( J) 

1  -0.108214600E  ♦  3 

2  -0.432858401E  ♦  3 

3  -0. 97393 1402E  ♦  3 

4  -0.173143360E  ♦  4 

5  -0.27053650  IE  +  4 


31  -0.103994231E  +  6  961.0000 

series  truncated  after  31  poles 


APPENDIX  B 


USE  OF  COMPLEX  ANALYSIS  FOR  CASE  OF  UNEQUAL  D,  AND  D2 

This  appendix  presents  the  poles  of  the  response  function  for  a 
double-layer  system  with  the  given  parameters  (see  Figs.  Ill -5  and  III- 
6).  Compare  to  Appendix  A. 

Case  I  (Table  III-1):  Dj  =  1.0  cm^/Msec 

D2  =  0.1 

H  =  1000  Msec-1 

_J  SN(J)  Multiplier 

1  -0.75911 1699E  2  1.0000 

2  -0. 343936 153E  ♦  3  4.5308 

3  -0. 84098 1695E  +  3  11.0785 

4  -0.155907427E  ♦  4  20.5381 

5  -0.233610574E  ♦  4  30.7742 


31  -0. 7583402 16E  +  5  998.9837 

series  truncated  after  35  poles 


APPENDIX  C 


PARTICLE-SIZE  ANALYSIS 

The  Le  Mont  Image  analyzer  was  used  to  obtain  the  following  particle- 
size  information. 

SAMPLE  ID:  SrCrO,, 

Number  Count;  148 


Measurements  for  Each  Physical  Parameter 


Area 

(ym2) 

Ratio  of 
Width 
to  length 

Width 

(yo) 

Length 

(ym) 

Angle 

(deg) 

Perimeter 

(ym) 

Avg 

4.39 

0.163 

0.7162 

5.2099 

82 

11.7 

RMS  Dev 

4.66 

0.075 

0.3448 

3.1623 

54 

6.79 

%  RMS  Dev 

106 

46 

48 

61 

66 

58 

Median 

2.97 

0.149 

0.6907 

4.2300 

79 

- 

Most 

Probable 

3.40 

0.125 

0.8150 

3.2500 

23 

- 

SAMPLE  ID: 

Ma0*Si0o 

Number  Count:  105 

Measurements 

for  Each  Physical  Parameter 

Area 
(yn  ) 

Ratio  of 
Width 
to  length 

Width 

(ym) 

Length 

(ym) 

Angle 

(deg) 

Perimeter 

(ym) 

Avg 

7.25 

0.262 

1.2728 

4.6207 

86 

11.4 

RMS  Dev 

9.17 

0.109 

0.9880 

2.5176 

54 

6.32 

%  RMS  Dev 

126 

42 

78 

54 

62 

56 

Median 

3.98 

0.236 

0.9794 

4.1078 

84 

- 

Most 

Probable 

7.30 

0.225 

0.8150 

5.1500 

164 

- 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for 
national  security  projects,  specializing  In  advanced  nllltary  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  Investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systeas.  Vital  to  the  success  of 
these  investigations  la  the  technical  staff's  wide-ranging  expertise  and  Its 
ability  to  stay  current  with  new  developaents.  This  expertise  Is  enhanced  by 
«  research  prograa  aiaed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systeas.  Contributing  their  capabilities  to  the 
research  effort  are  these  Individual  laboratories: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynaalcs;  chemical  and  electric  propulsion,  propellant 
chealatry,  chemical  dynaalcs,  envlronaental  chealstry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  atructural 
control;  high  temperature  thermomechanics ,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exclmer  laser  development  including  chemical  kinetics, 
spectroscopy,  optlcsl  resonators,  beasi  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory :  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of -fleld-of -view  rejection, 
applied  laser  spectroscopy,  laser  chealstry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vecuua  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
envlronaental  chemistry. 

Computer  Science  Laboratory:  Prograa  verification,  prograa  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  Intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solld-stste  lasers,  optlcsl  propsgatlon  and  communications; 
microwave  semiconductor  devices,  microwave/ millimeter  wave  measurements, 
dlegnoettcs  and  radiometry,  alcrowave/nllllmeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  In  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospherlc,  auroral  and  cosmic  ray 
physics,  wave-particle  Interactions,  magnetospherlc  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  uaif^  atmospheric  radiation;  solar  physics.  Infrared  astronomy. 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
lnstrtaaentatlon. 


